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INTENSITIES, POLARIZATION AND ELECTRON DENSITY 
OF THE SOLAR CORONA FROM PHOTOGRAPHS TAKEN AT 
THE TOTAL SOLAR ECLIPSE OF 1952 FEBRUARY 25 


H. von Kliiber 
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Summary 

The corona was photographed at the total solar eclipse of 1952 February 25 
with a quadruple camera of focal length 200 cm. Six plates were obtained, 
each with four corona images, exposed simultaneously, three through pola- 
roids with planes of polarization 60° apart, and one without polaroid. Each 
plate carried relative and absolute calibrations. An extensive photometry 
of these plates has provided intensities, amount and direction of polarization, 
a separation of K and F coronas, and electron densities. The results are 
summarized in graphs. Comparison is made with van de Hulst’s model 
and with other observations made at the same eclipse. 





General remarks.—Recently van de Hulst (1) computed a general model 
of the solar corona, based mainly on the exhaustive studies of Baumbach (2), 
who used all available eclipse observations from 1go1 to 1932. Van de Hulst’s 
model is considered by most astronomers to be the best obtained so far. The 
purpose of the present work was to produce under favourable conditions a 
comprehensive set of measures of intensity and polarization of one particular 
corona, and thence to derive the electron density and the now generally accepted 
division into K and F components (W. Grotrian (3)), for comparison with the 
van de Hulst model and for use in future theoretical investigations. 

The eclipse of 1952 February 25, observable in the desert climate of the Sudan, 
offered an unusually good opportunity for photometric measurements of the kind 
required. Furthermore, since many other expeditions and observers, both in 
the Sudan and elsewhere, made measures at the same eclipse (4), the structure 
of the corona on that particular day is probably better known than on any other 
occasion so far. 

The instrument.—For photometric reasons it was considered essential to make 
intensity and polarization measures together in groups, each group strictly 
simultaneously. An arrangement proposed long ago by Fesenkov (5) was 
used, the corona being imaged through three different lenses, each provided 
with a polaroid screen, the planes of polarization being oriented 120° apart. 
A fourth lens photographed the corona in unpolarized light. The method allows 
both the amount and the direction of polarization to be obtained independently. 

A specially constructed quadruple camera was used, with four identical 
lenses (cemented doublets) of f=200 cm, 6=7°5 cm, arranged at the four corners 
of a square, the diagonal separation of their centres being 15cm. ‘The lenses 
imaged simultaneously on one 10 in. x 12 in. photographic plate. The useful 
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field was a circle of radius at least 4 solar radii from the centre of each image, in 
some directions several radii more. It is important to have available as much 
field as possible to measure the sky fog at sufficiently large distances from the 
Sun. The fourth lens, working without polaroid, was stopped down so that its 
image was of about the same intensity as those photographed through polaroid. 

No polaroid was available of sufficiently good optical quality to be used over 
the lenses. After a large number of tests, the polaroids were therefore mounted 
about 20 cm in front of the plate, a position giving a satisfactory compromise 
between the best optical definition and the least scattered light. Every possible 
precaution was taken to avoid stray light; the camera was divided into four so 
as to isolate each light beam, stops were provided at various places, and the 
whole interior was painted with a special matt black paint. 

The most desirable way of using such a camera would be on an equatorial 
mounting, but this was impossible for technical reasons. Instead it was mounted 
horizontally and fed by a 16-in. coelostat with an aluminized mirror. This 
procedure may be viewed with some doubt because it is well known that reflection 
at a metal surface alters the amount and character of the polarization. <A careful 
study, by both calculation and experiment, was therefore made to estimate the 
size of any effect on our results. The angle of incidence on the coelostat mirror 
at the moment of totality was 38°. The degree and form of the polarization 
introduced were calculated for this special case with the usual formulae (6), 
using the optical constants for aluminium published elsewhere (7). The result 
was that with completely unpolarized incident light, the reflected beam would 
show very slight elliptical polarization, equivalent to a plane polarization of just 
under 2 percent. This value agrees with recent calculations by J. N. Hodgson (8), 
who also gives useful approximation formulae for such work. 

Furthermore, this result was confirmed by direct observation with the actual 
mirror as it was used during the eclipse. A large piece of cardboard was 
illuminated from a distant electric lamp in such a way that no fringes could be 
seen when carefully examined by a Savart prism. The same cardboard was 
then viewed by reflection from the coelostat mirror, again using the Savart prism. 
No fringes could be seen with certainty for angles of incidence less than 40°. 
On the other hand, Savart fringes were seen when viewing the cardboard directly 
through a glass plate inclined about 25° from the normal and introducing about 
2 per cent polarization (g). This means that polarization greater than 2 per cent 
could be detected with certainty by the Savart prism (10), and that the amount 
introduced by the coelostat mirror at the eclipse was no greater than this. Since 
the present photographic investigation was not expected to give errors smaller 
than 2 per cent, the polarization introduced by the mirror has been neglected. 
The labour involved in removing the small mirror polarization from many 
thousands of measurements would have been quite out of proportion to the 
negligible improvement following in the results. 

Since the degree of polarization and its orientation are to be determined by 
combining intensity measurements from three cameras equipped with polarizers 
it is extremely important that the photometric relations of the cameras should 
be known accurately. Workers in this problem have in the past usually assumed 
the corona to be plane-polarized radially, so that the intensity as viewed by any 
two cameras should be equal at any point along a radius dividing the angle between 
their two polaroids. But this procedure could not be applied here, for we wished 
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to make no assumption about the direction of polarization. Thus special tests 
were carried out, partly before and after the eclipse in the laboratory, and partly 
at the eclipse site itself, to find the corrections required to reduce each camera 
toacommon scale. This was done first by photographing a uniformly illuminated 
sheet of cardboard, sufficiently free from polarization, using the cameras both 
with and without their polaroids. Also, the Moon was photographed with 
adjustments close to those used for the eclipse. A completely overcast sky was 
photographed twice, choosing a part on which no Savart fringes could be seen. 
Further, on the night following the eclipse, star fields and also trails of bright 
stars were photographed with all cameras. 

All these methods showed good agreement with each other. Only very small 
relative corrections were indicated and the size of the corrections which were 
actually applied to all the measurements appear to be well established. The 
transmission of the polaroids, tested separately, was found to be sufficiently 
uniform for present purposes over the whole field. In the wave-length region 
used, a pair of crossed polaroids transmitted about 1°5 per cent. 

The observations.—By courtesy of the military authorities, the observations 
were made at Fort Stanley, Khartoum. The instrument was set up in an open 
courtyard within the Fort, where it was well protected from wind and sand. Water 
and electricity were available, and a sufficiently large darkroom could be improvised 
nearby as well, an important point for photometry, as will be explained below. 

The local circumstances of the eclipse were: 


TABLE I 
Geographical coordinates of Fort Stanley A=32° 32’6 E 
= 15° 363 N 
Height above sea level 380m 
2. contact 9" om 158-5 U.T. 
3. contact gh 12M 248-5 U.T. 
Duration of totality 1898 
Position angle of 2. contact 39° 
Position angle of 3. contact 233 
Height of Sun at mid-totality * 62° 
Sun’s apparent diameter a3 
Moon’s apparent diameter 33°°5 
Position angle of the Sun’s axis —20°°3 
Hel. latitude of the centre of the Sun’s disk —7'1 


Table II gives details of the exposures. These were made with different 
stops and exposure times, arranged so that the whole corona from the limb 
to about 4 radii from the centre was covered with images of a favourable range of 
density for photometry and with sufficient overlap. The lens stops were on 
one single metal plate and needed to be withdrawn once only during totality, 
a change requiring less than one second. 

Figs. 1 and 2 are reproductions from our plates 2 and 4; many faint details 
are lost, as usual, in the reproductions. 

A good deal of thought had been given to selection of a suitable photographic 
emulsion, following laboratory experiments. A moderately high y was wanted 
and, in view of the large intensity range in the corona, as much protection as 
possible against halation. This is believed to be very important in this kind of 
work, although it seems often to be neglected, and more attention should be paid 
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to it at future eclipses. ‘The Kodak Research Laboratories, Harrow, kindly 
prepared several test emulsions for use, using various anti-halation dyes. These 
emulsions were tested at the Cambridge Observatories, using a long range step- 
wedge and up to as much as 10™ over-exposure. The problem was to decide 
on the dye concentration, too much dye causing too much loss of sensitivity. 
Eventually Kodak experimental emulsion 50707 SDC 2588, tropical-hardened, 
was selected. This is a blue sensitive I-O emulsion to which a yellow dye has 
been added. Fig. 2, taken from the long exposure 4, shows that, even with a heavy 
over-exposure for the innermost corona, very little halation is visible. 


TABLE II 
Exposed Exposure Apertures Aperture 
Plate after sec. for the of the Emulsion 
2. contact 3 lenses __ free lens 
sec. with cm 
polaroids 
cm 

I 5-9 4 2°5 1°4 Kodak 
2 30—46 16 2°5 1*4 Kodak 
3 60 —96 30 7S 4°2 Ilford 
4 115 —145 30 7's 4°2 Kodak 
5 162 — 166 4 o-% 4°2 Kodak 
6 178 — 184 6 7€ 4°2 Kodak 


The colour sensitivity is that given by Kodak (11) for an O-type emulsion, 
extending from the violet to about 4900 A. Our photometry thus refers to a 
fairly wide wave-length range with a centre of gravity at about 4300 A. 

For exposure No. 3 an Ilford tropical plate (batch S 2466), having about the 
same colour sensitivity and also anti-halation protection, was used. ‘This plate, 
about five times faster than the Kodak plates, was chosen to show the outer corona 
and to obtain values for the general sky light, in case exposures 2 and 4 should 
failto doso. Exposure 6 was a reserve exposure, fitted in near the end of totality, 
and actually has not been needed for the complete series. 

The sky was remarkably good during totality, conditions being very favourable 
for photometric work, and the programme was carried through without any 
hitch. The change of stops and the exposures were carried out by hand by 
Mrs L. von Kliiber. Many trials had previously shown that all four cameras 
received exposures equal to within o-1 sec and that no systematic error was to 
be feared on this point. 

Calibration.—We believe it necessary that in this kind of work each plate 
should carry its own intensity marks. As in earlier investigations of the author 
(12) the calibrations were carried out quite simply but in a manner which 
is believed to give high relative accuracy over the whole corona and at the same 
time absolute intensities good to 10 or 20 percent. This figure must be regarded 
as satisfactory when one recalls the many sources of error inevitably involved 
in all photometry of this nature. 

With calibration needs in mind the instrument was set up in Fort Stanley 
only a few yards from a disused storeroom which could easily be converted into 
a good and entirely light-proof darkroom. Sunlight was reflected from the 
coelostat through a hole in the wall directly into the darkroom. The diameter 
of the hole could be varied at will by means of calibrated circular metal stops. 
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Fic. 2.—Corona as taken with the quadruple corona camera at the eclipse of 1952 February 25 
at Khartoum. 


f=200 cm o=420cm exp. 308 


H. von Kliiber, Intensities, polarization and electron density of the solar corona. 
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For calibration each photographic plate was transferred to a specially prepared 
auxiliary plateholder into which had been fitted four step-filters (provided by 
the Kodak Research Laboratories) about 3 cm from the edge of the plate, each 
adjacent to one image of the corona. The step-filters each had twenty steps, 
3mm high, the interval between two successive steps being o-1 in density 
(=o0™-25), so that in all an intensity range of 1: 100 was covered. The sunlight 
entering the darkroom was reflected by a diffusing surface of chalk and at a fixed 
known distance from this the plate was exposed through the step filters, the 
exposure time in each case being exactly the same as for the corona. In the case 
of corona images obtained through polaroid the calibration was also made through 
polaroid. One may safely assume that the spectral distribution of the light of 
the Sun and the corona for this work are sufficiently nearly identical, so that the 
calibrations were made with both the correct colour distribution and the correct 
exposure time. 

Further, taking into account the lens aperture, the loss of light at the lenses 
when used for the corona photographs, the diameter of the entrance aperture 
in the darkroom wall, the reflectivity of the chalk (assumed to be g5 per cent (13)), 
and the distance between chalk and plateholder (known to 1 per cent), an 
elementary calculation (14) gave the light intensity incident on each step-filter 
in units of the mean brightness of the Sun’s surface. This unit has been used 
for all the corona intensity measurements in this paper. Care was taken, of 
course, to print the calibration marks with the Sun at the same zenith distance 
as it was during totality. In addition to some check calibrations taken immediately 
after the eclipse, a large number of test calibrations made on other days, both 
before and after the eclipse, showed remarkably good internal consistency. 
Provided days with well-marked dust and haze were avoided, the sky in the desert 
climate of the Sudan gave very consistent transparency. 

The corona exposures were arranged in such a way that the photometrically 
best ranges overlapped, and then the measures from each could be linked together 
to give one consistent system of relative intensities over the whole corona. 
Furthermore, the calibration just outlined allowed the zero point of the intensity 
scale for each plate to be fixed independently in units of the average surface 
brightness of the Sun’s disk. Had the calibrations been faultless these zero points 
would have been confirmed exactly when fitting together the measures from each 
plate. For the four plates ultimately used the following corrections were 





required : 
Plate 1 +0-04 in log / 
2 — 0°02 
4 — 0°04 
5 +0°02 


Considering the large relative uncertainties to be expected in such photometry, 
this must be regarded as very good agreement. ‘The photometric system 
established in this work should be regarded as observationally quite reliable. 

The photometric scale of course depends critically on the calibration of the 
step-filters themselves. Great care was taken with this. Most of the calibrations 
actually used were made before the eclipse in the laboratory and have been 
repeated since. These calibrations were carried out photographically in three 
independent ways, using the same plate emulsion as for the eclipse: com- 
parisons made were (a) with a tube sensitometer, (b) with a Hilger rhodium 
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step-filter, calibrated by the N.P.L., and (c) with the inverse square law. The 
three methods gave very good interagreement. 

Development.—The plates were brush developed (15) in Khartoum to a fairly 
high gamma, using undiluted Kodak D1gb developer for five minutes (16). 
This was followed by a stop bath (Kodak SB4) and Kodak acid fixer with 
hardener. Although the daytime temperature outside on most days went to 
35 °C or more, no difficulty in processing was encountered. In the very dry 
Sudan air the processing solutions could easily be kept down to 18 °C by storing 
the bottles in water-filled unglazed earthenware pots of the type widely used 
by the Sudanese for cooling liquids. Drying the plates needed no special 
precautions. All plates were found to be of excellent quality for photometric 
work. 

Photometry.—To keep the very large amount of photometric work within 
reasonable limits only plates 1, 2, 4 and 5 were measured. ‘The spare plate 6 
would have contributed little more. Plate 4 showed just about the optimum 
density for measurement of the background of sky and scattered light. Plate 3 
was unnecessarily dense for this purpose. 

The correct orientation of a corona picture with regard to position angle often 
meets with some difficulty. Here three methods were used: 

(a) On the night after the eclipse, with the camera in the same position as 
during totality, star trails were photographed on an auxiliary plate. 

(b) The longer eclipse exposures showed several fine star images, especially 
BD 5850, 5855 and 5888, which not only proved that guiding had been excellent 
during totality, but also provided a very good orientation. 

(c) All plates showed several prominences at the Sun’s limb, the two most 
conspicuous being indicated in Fig. 5. Their position angles could be taken 
from observations made at about the same time at the Wendelstein, Ziirich and 
Locarno Observatories. Several different determinations showed that the 
position angle of the prominence near the W limb was so near 270° that for this 
work it was sufficient to assume it to be exactly at 270° with respect to the Sun’s 
axis (17). 

Unfortunately no photometer was available which would give direct traces 
of isophotes or of densities in concentric circles around the Moon, as would have 
been convenient for this work. Instead, the Cambridge Observatories’ micro- 
photometer was used. Its slit dimensions, as projected on the plate, were 
0°15 x 0-02 mm, (0’-22 x 0’-03 on the sky), the long side being tangential to the 
Moon’s limb. Straight sections through the centre of the Moon were recorded 
at all position angles in steps of 5° or 10°. 

The main uncertainties in our final values are believed to be due to this 
photometry in radial strips. To deduce the polarization, three measures from 
three different images had to be combined. It is easy to see that even a very 
small error in the assignment of the position angles would introduce small local 
differences between the three measures. This is especially true of radial streamers 
for which the density gradient is high and for which this measurement procedure 
is unsuitable. For this reason we have intentionally omitted some values around 
the strong streamer at about position angle 80°. It is hoped to make a special 
study of this streamer later. Frequent repetitions of some of the measurements 
showed good internal agreement, which suggested that the radial tracing method 
is sufficiently accurate for the present purpose of getting a good general picture 
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of the corona as a whole. Making the microphotometer records required many 
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weeks. The intensity steps printed on the plates were measured frequently, 
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sky fog found at about p=6. The other less exposed plates do not reach the sky 
fog and therefore show nothing of this kind. This effect can hardly be explained 
simply by scattering of light from the bright corona ring or by photographic 
effects because it is so uniform over the whole disk of the Moon. Actually we 
have very little accurate information about the earth-shine on the Moon at new 
moon. The best values are perhaps those extrapolated from observations of 
Danjon (21) from which for the new Moon in our units an average value of the 
earth-shine of logo-4 to logo-5 can be derived, so that we should expect to find 
the intensity of sky fog plus earth-shine in Fig. 4 to about log1-3 instead of 
log 1°5. On the other hand it is well known that the albedo of the Earth depends 
mainly upon fortuitous cloud conditions over large areas and varies appreciably 
(at least between 0-32 and 0-52 according to (21), but probably more according 
to (22)). Obviously the earth-shine must vary in a similar way. Unfortunately, 
very little and not particularly accurate quantitative information seems to be 
available of the earth-shine from Sun eclipse observations (23, 24,25). Since 
the sky light on our plate was found to be already fairly weak at p = 6, the additional 
intensity on the Moon’s disk could be explained only by the not very plausible 
assumption of a fairly faint halo, quite uniform over the Moon’s disk and then 
decreasing outwards fairly rapidly. If one tried to introduce a correction for this 
hypothetical effect, the values for log (K+ F), as given below in Figs. 11 and 12, 
would not be affected up to a distance of p<2-2; from here on small negative 
corrections would probably be needed which near the end of our range at p ~ 3°5 
might attain —o-1 to —o-2 in log(K+¥F). The same correction would alter 
our values for polarization only by a second-order effect in the outermost parts 
of the corona. But until the cause for the unexpectedly high intensity on the 
Moon’s image (also noticed by other observers) is properly understood we have 
thought it better not to introduce a dubious correction for the effect. 

All our measures were corrected for scattered light by a method developed 
long ago by K. Schwarzschild (26), assuming that the sky light is scattered 
uniformly over the fairly small field (2° x 2°) covered by our instrument. This 
seems plausible, for the chief contribution to the sky light comes by scattering 
from outside the belt of totality. A number of laboratory tests showed that the 
instrumental scattering from coelostat mirror etc. was sufficiently uniform. 

Throughout this paper we shall use the following notation: K=intensity of 
the electron corona; F=intensity of the dust corona; (K+ F)=observed 
intensity of the combined Kand F coronas; ~:, , », = directly observed polarization 
of the K+¥F coronas’; p,=polarization of the K corona. 

All intensities are given in units 10~!° of the average brightness of the Sun’s 
disk. 

Position angles are referred to the rotation axis of the Sun, measuring from 
north through east. 

Positions within the corona are given in solar radii from the Sun’s centre; 
r is the distance within the real three-dimensional corona, p is the distance in 
the projection seen on the sky. 

In many investigations the corona has been divided into polar and equatorial 
regions respectively. ‘To achieve the main aim of this paper, namely a direct 
comparison with van de Hulst’s figures (27), the ‘‘ polar corona’’ was taken 
from results averaged over position angles 330° to 30° and 150° to 210°, while 
the ‘‘equatorial corona’’ was taken from position angles 40° to 140° and 
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Fic. 5.—Isophotes in log (K+-F) as observed; inner region of the corona. The positions of two 
prominences are indicated by P. 








Fic. 6.—Isophotes in log (K+ F) as observed; outer region of the corona. 
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220° to 320°. This crude dichotomy should not be regarded as always adequate 
for discussion of the corona problem. 

After correcting for scattered light, the intensities from the photographs 
taken without polaroid could be combined to give the distribution of surface 
brightness of the K+F corona, expressed in absolute units. The results are 
shown in Figs. 5 and 6. In the case of the photographs taken through polaroid, 
the correction to be applied was for the appropriate component of the linearly 
polarized sky light. In the way described below, three sky measurements were 
combined to give the amount and direction of the sky polarization. Near the 
Sun this was found to be 11 per cent, inclined at 9° to the vertical. This angle is 
not regarded as very reliable, because it is determined as the resultant of three 
rather small vectors and a small error in one of these could alter the angle quite 
appreciably. The figures are nevertheless in good agreement with those of 
Michard (19) from the same eclipse and the same observing site, who found 
a polarization 12 to 14 per cent oriented approximately vertically. Our values are 
also in good agreement with those of Ohman (28), obtained in 1945 using the 
same method. 

After applying the necessary sky corrections there were three separate measures 
for each point in the corona which had to be combined to give the amount and 
direction of the polarization. By definition the amount of polarization is 

Ba 0) 
where J, and J, are the intensities of the tangential and radial components (electric 
vector) respectively. Denoting the position angle of the plane of vibration by « 


and the three measured intensities by a, b and c, we get for polarizers orientated 
120° apart (28), 





a=TI,cos*«+TI,sin*«, 

b=I,cos* (120° —a)+TJ,sin? (120° —«), 

c=I,cos? (60° —«) + I, sin? (60° — «). (2) 
The amount of polarization is given by 


2\/a(a—b)+b(b—c)+c(c—a) 











Petr) = aahee (3) 
and the direction relative to the a component by 
—b 
tana= 73 —— —: (4) 


The plane of vibration of one of the three polaroids was adjusted nearly 
parallel to the Sun’s axis, and was determined within 1°, using a special small 
device with an auxiliary polarizer and a plumb-line. 

The equations given above had to be solved for a considerable number of 
points along each of our radial records. The results were plotted graphically, 
allowing all amounts and directions of polarization to be read off easily. ‘The 
large number of detailed figures need not be published here; they are suitably 
condensed in Figs. 7 and 8, from which it is hoped that all desired information 
can be extracted. We have deliberately avoided giving values around position 
angles 70° owing to the very high polarization in the bright streamer there. 
The present procedure of radial measurement is not entirely satisfactory when 
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Fic. 8.—Direction of polarization (magnetic light vector) as observed. The length of the vectors 
indicate average relative intensities of polarization. (See also Fig. 13.) 
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dealing with such a feature. Fig. 7 gives contours of equal polarization over the 
corona. The directions of the magnetic vectors are plotted in Fig. 8, which gives 
the mean values for concentric circles at p= 1-2, p=1-6 and p=2°4. The lengths 
of the vectors at the same time give the average relative intensities of polari- 
zation as found at these distances (see also Fig. 13). 

There is a scatter in the directions of these vectors, rising in a few places to 
as much as 10° from the radial direction. This is not surprising, since each vector 
here is the resultant of three parent vectors. On the whole the direction of 
polarization is practically radial and, apart perhaps from a small irregularity in 
the NW quadrant, possibly connected with the rather high local polarization 
found there, it is doubtful whether the deviations from the radial directions are 
significant. 

Determination of electron density.—Having now a large number of values of 
the intensities and polarizations, which so far all refer to the sum of the K and F 
corona, it is possible to deduce electron densities for the K corona. We must 
keep in mind that the procedure now to be followed assumes a radially symmetrical 
corona. We know that this is not strictly the case and deviations from this 
assumption must be allowed for, especially near the poles and long streamers. 
Without more or less complicated and uncertain assumptions we can hope only 
to get fairly good average values for this particular corona. However, these 
may not be so far from the real values; it is well known that despite the apparently 
complicated structure of the corona on the photographic plate (or still more 
when viewed by the naked eye) a careful photometry usually gives remarkably 
smooth isophotes. This may be taken as an indication that average values are 
not likely to be much affected by local structure. 

To derive electron densities with the generally accepted assumption that the 
light of the K component originates from scattering by free electrons, a rather 
complicated integration is necessary. Here we have used a method developed 
long ago by K. Schwarzschild (29), and used later by Baumbach (2) but in the 
form employed by van de Hulst (1), to whose work reference should be made 
for details of the mathematical treatment. The procedure takes the polarization 
into account, and thus the fact that the observed intensities are the sum of the 
polarized K and the probably non-polarized F component. 

Using still the indices t and r respectively for the tangential and radial intensity 
components, we have by definition 


| (K+FP). dg. m=Pe=Ki-K, ( 

where ~,,, ») and (K+F) are known from our observations. 
On the other hand, if K,— K, can be expressed by a power series of the form 
K,-—K,= > kyr-, (6) 


wn 
~ 


where r is the distance of the point under consideration from the centre of the 
Sun, the electron density N(r) can be represented by another power series 
of the form 


F k, : I 
N(r)= .. , ——. —. .—— , " 
= 20 FC AG)-Be) ” 
Here the coefficients a, can be obtained in terms of gamma functions (30) ; 
they are also given by a graph in (25) or in a table in (1). The constant C takes 


into account the scattering of light from free electrons. If the intensities are to 
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be in units of 107° x average brightness of the Sun, C=3°44x10-4cm. The 
coefficients A and B take into account the darkening towards the limb of the Sun 
and are also given in a table in (1). Following van de Hulst we have accepted 
A and B values corresponding to a coefficient of limb darkening 0°75, which 
corresponds near enough to the effective wave-length of our photometry. Using 
formula (7) electron densities could be calculated for each of our sections across 
the corona. Again to avoid much tabulation we give the results graphically, in 
Figs. g and 10, which show contours of equai electron density. 

Similarly the values for the intensity of the K component alone could be 
found. Following (1) we have the relation 





K,= > hk, .r-* (8) 
and 
‘ h, 1 
N(r). A(r).C= Y= .—. (9) 
37 as—y 
Having now the values for N(r) a power series can again be found of the form 
: h, 
N(r). A(r).C= » (10) 
and thence we find immediately 
. h, 
Km Sha (11) 


With both K,—K, and K, the values for K itself are known. After very 
slight smoothing they are represented further below graphically in Figs. 16 and 17. 

With K determined and with (K+ F) from the observations, both F and p, 
could immediately be found for all positions at which measures were available. 
As was to be expected, the values for log F show no marked systematic dependence 
on position angle. Thus they could be averaged and represented as a function 
of distance only. There should be no dependence of F on the phase of the solar 
cycle, so that the results have been compared with the corresponding van de Hulst 
figures, shown further below in Fig. 19. It is estimated that our values could 
easily be in error by 10 to 20 per cent, so that the agreement is really remarkably 
good; there is a tendency for our values to be slightly lower than van de Hulst’s, 
but the reality of this difference is questionable. 

If a few abnormal regions with strong streamers are excluded, the polarization 
px is fairly uniformly distributed in position angle. Further below Fig. 18 gives 
a comparison with the van de Hulst model. There is a tendency for the measured 
polarization to be lower than that of the model beyond 2 radii from the centre. 
In comparing the measured and computed values it should be remembered that 
measurement errors inevitably tend to be greater as one goes outwards from the 
Sun. 

In a symmetrical K corona in which the decrease of intensity outwards follows 
the law 
< (12) 
p, would be expected at large distances to tend ton/(n+2). ‘Taking m from the 
intensity measures one would then get a limiting polarization of about 70 per cent. 
Following similar considerations Michard (19g) obtained this value also for the 
same eclipse. The measured values tend to about 60 per cent. 


K= 
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Fic. 9.—Lines of equal electron density in log N calculated from observations; inner region 
of the corona, 





Fic. 10.—Lines of equal electron density in log N calculated from observations: outer region 
of the corona, 
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One of the chief aims of this investigation is to compare our results with the 
values of intensity and polarization given by the van de Hulst model. This is 
done in Figs. 11-19, polar and equatorial regions being considered separately. 
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Fic. 11.—Log of intensities of (K+-F) against solar radii, for polar region: 
(1) Cambridge observations 
(2) Model van de Hulst, minimum (27) 
(3) Model van de Hulst, maximum (27) 
(4) N.R.L. observations (33) 
(5) Michelson’s observations (47) 
(6) Michard’s observations (31). 
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Fic. 12.—Log of intensities of the (K+-F) component against solar radii, for equatorial region: 
(1) Cambridge observations. 
(2) Model van de Hulst, minimum (27) 
(3) Model van de Hulst, maximum (27) 
(4) N.R.L. observation (33) 
(5) Michelson’s observations (47) 
(6) Michard’s observations (31) 
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Figs. 11 and 12 show that our values for (K + F) are in fairly good agreement 
with the model. The slight tendency for the measured values to be higher very 
near to the Sun’s limb must not be regarded as necessarily real for two reasons: 
(a) the gradient there is very steep and cannot be determined very accurately 
on fairly small-scale photographs such as have been used here, and (d) parts of 
the innermost corona were behind the Moon’s disk when our photographs were 
taken. From about 2:5 radii outwards our values tend to lie systematically 
below those of the model. This cannot be adequately explained by a possible 
over-correction for scattered light. 

Fig. 13 shows that the measured polarization in the equatorial region follows 
the ‘‘minimum’”’ model corona well. In the polar region the measures lie well 
between the curves for ‘‘maximum’”’ and ‘‘minimum’”’ model coronas. In 
both regions the maximum polarization is found near r= 1-22. 
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Fic. 13.—Polarization of the (K-+-F) component against solar radii: 
(1) Cambridge observations, equatorial region 
(2) Cambridge observations, polar region. 
(3) Model van de Hulst, maximum for pole and equator (27) 
(4) Model van de Hulst, minimum equator (27) 
(5) Model van de Hulst, minimum pole (27). 


One of the most interesting physical quantities connected with the corona, 
the electron density, is found for the polar region (Fig. 14) to have values well 
between those for ‘‘maximum”’ and ‘‘minimum”’ models. In the equatorial 
regions (Fig. 15), for r>2-2, our values are slightly greater than those of the 
model and also show a somewhat steeper gradient than in the model. Again 
the values very near the solar limb are somewhat uncertain, for the reasons 
mentioned above. 

Taking the K corona alone (Figs. 16 and 17) the polar region measures 
again lie well between the ‘‘maximum”’ and ‘‘minimum”’ models; for the 
equator they agree well with the ‘“‘ maximum ’”’ model. 

The polarization calculated for the K component alone (Fig. 18) has been 
discussed above. As there were hardly any differences between pole and equator 
only one curve is shown here. The limiting value for large r is only a few per cent 
lower than that for the model. 

Our figures for the F corona, mentioned already above and represented here 
in Fig. 19, agree very well with van de Hulst’s values. 

Remembering that the model corona is necessarily a somewhat ideal picture, 
we may say that our 1952 measures on the whole agree well with the model, 
with very little certain systematic differences. 

Comparison with other observations of the same corona.—As we have expected 
and hoped, many corona observations from this eclipse have been published, 
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Fic. 14.—Log of electron density N of the K component against solar radii, for polar region: 
(1) Computed from Cambridge observations. 
(2) Model van de Hulst, minimum (27) 
(3) Model van de Hulst, maximum (27) 
(4) Michard calculated (19, 31). 
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Fic. 15.—Log of electron density N of the K component against solar radii, for the equatorial 
region: 
(1) Computed from Cambridge observations 
(2) Model van de Hulst, minimum (27) 
(3) Model van de Hulst, maximum (27) 
(4) Michard calculated (19, 31). 
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Fic, 16.—Log of intensity of the K component against solar radti, for the polar region: 
(1) Computed from Cambridge observations 
(2) Model van de Hulst, minimum (27) 
(3) Model van de Hulst, maximum (27). 
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Fic. 17.—Log of intensity of the K component against solar radii, for the equatorial region: 
(1) Computed from Cambridge observations 
(2) Model van de Hulst, minimum (27) 
(3) Model van de Hulst, maximum (27). 
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Fic. 18.—Polarization of the K component against solar radii: 
(1) Computed from Cambridge observations 
(2) Model van de Hulst, minimum pole (27) 
(3) Model van de Hulst, minimum equator and maximum pole and equator (27). 
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Fic. 19.—Log of intensity of the F component against solar radii: 
(1) Computed from Cambridge observations 
(2) Model van de Hulst (27). 


so that taking everything together we have rather good data concerning the 
corona on this particular occasion. Among the measurements which have 
appeared so far are a number dealing with the surface brightness and polarization 
of the corona which may be compared briefly with our observations. 


(a) The French expedition at Khartoum (19, 31) 


Michard et al. carried out polarimetric and absolute photometric measure- 
ments during this eclipse. Comparison with the present work is restricted because 
they worked in the red region of the spectrum and covered mainly the range 
2°5 <p <20, so that we have only a few overlapping values. The few figures for 
polarization that we have in common are given underlined in Fig. 7 and show 
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good agreement. Reference was made earlier to Michard’s measures of sky 
brightness and polarization. His intensity measurements for (K+) provide 
a valuable extension to those of this paper (see Figs. 11 and 12). The gradients 
of the French measurements are somewhat flatter than ours, at both equator 
and pole. This difference might be associated with the F corona following the 
different wave-length region used but it persists in the distribution of electron 
density deduced for the K corona only (see Figs. 14 and 15). The French 
measurements have in turn been extended (at the 1954 eclipse) up to over 50 p 
by Blackwell (32). 


(b) The United States Naval Research observations at Khartoum (33) 

These measurements were a relative photometry only; they included no 
absolute calibration and no polarization. They were made with the Gardner 
camera of 6m focal length and cover a range 1<p<2:5 in the corona. The 
photographs were taken on Kodak Panatomic X emulsion and therefore 
presumably cover the whole visible spectrum to about 6300A. Since the 
plates were calibrated only for relative intensities, only their slope can be compared 
with our values. The absolute values given in (33) were obtained by adjustment 
to van de Hulst’s model, which was also used for separation of the K corona 
since no original measurements for the polarization were available. 

If instead of using the van de Hulst model the N.R.L. figures are adjusted 
to our absolute values, a shift of about 0-14 in log (Ff + K) is required and there 
is then very good agreement with our curves (see Figs. 11 and 12), which in view 
of the considerable differences in telescopes, calibration and photometric 
procedure may be considered very satisfactory. 


{c) Italian expedition at Khartoum (34, 35, 36) 

The Italian expedition from Arcetri, also at Khartoum, obtained a set of 
photographs through polaroids, using a double camera (f=250 cm), covering 
the range 1:08<p<i-5. The resulting measurements have been derived from 
different photographs with different effective wave-lengths and are not strictly 
comparable with ours. Michard (31) has reduced them to a more or less 
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Fic. 20.—Polarization of the (K+F) component against solar radii, as observed by other 

expeditions at the same eclipse: 

(1) Cambridge observations, equatorial region 

(2) Cambridge observations, polar region 

(3) French and Italian observations, equatorial region (31) 

(4) French and Italian observations, polar region (31) 

(5) Vashakydze observations, equatorial region (46) 

(6) Vashakydze observations, polar region (46). 
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homogeneous system with his own measures. These values are to be found 
y 8 Apter it 
e in Fig. 20; they indicate a somewhat lower polarization than ours in the region 
: I'2<p<2°4. 
7 A further photometric and polarimetric study of the corona was made at 
e Khartoum by M. Cimino, Monto Mario (37), using a triplet camera (f= 30 cm). 
n No polarization values appear to have been published so far, but the mean values 
h for (K+ F) for the whole corona, shown in Fig. 3, reference (375), have been 
p converted to our scale here by bringing the values for p=2-o arbitrarily into 
coincidence and are shown in Fig. 21. As can be seen immediately for p< 2-0 
the Cimono values deviate strongly from ours in the sense of having a far smaller 
gradient, so that at p=3:0 the difference in log (K+ F) already exceeds 0°5. 
0 
rt (d) Russian expeditions 
e | Several Russian expeditions working in the eastern part of the track made 
e photometric and polarimetric measures of the corona at the same eclipse (38-47). 
e 
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Fic. 21.—Log of intensity of the (K+-F) component against solar radii; means over the 
| whole corona. 

(1) Cambridge observations 

(2) Cimono’s observations (37) 

(3) Nikolsky’s observations (44). 


Among the resulting publications is a very interesting paper by M. A. Vashakydze 
| (46) describing some work by the Abastumani Solar Eclipse Expedition to Cheeli. 
| He gives a careful polarimetric study of the corona carried out by a method very 
similar to that of the present investigation. He used a triple camera (f= 100cm), 
with three polaroids oriented 120° apart, and worked in the blue-violet. The 
calibrations (relative only) were made by a tube sensitometer and his four plates, 
each with three images and with their ranges overlapping well, covered a total 
range about 1:1<p<2°4. While the general picture of the distribution of 
polarization over the corona is in quite good agreement with ours, single values 
often show a fair amount of difference. ‘To compare the investigations as a whole 
| we have plotted the mean values of polarization for the polar and equatorial 
regions (Vashakydze’s Table 2) in Fig. 20. These values were averaged in 
position angle in the same way as ours and are therefore directly comparable. 


_aeewatre mes 





XUM 





222 H. von Kliiber, Intensities, polarization Vol. 118 


The curves are similar to ours in form, but in the region around p=1-4 there 
are systematic differences, our values being several per cent higher. Vashakydze 
has also calculated the directions of the polarization vector all round the Sun 
and the angles he finds are in excellent agreement with our values as shown 
in our Fig. 8. 

Another photometric investigation of the inner corona, with a focal length 
of 505 cm, and apparently on ordinary blue-sensitive plates, was carried out 
during the same eclipse by M. A. Michelson (47) at a Russian station near Cheeli. 
Unfortunately, difficulties were encountered in the absolute calibration and only 
relative values are available. Using the data in Michelson’s Table 5, and adjusting 
his relative scale to fit our absolute values as well as possible, we get the comparison 
shown in our Figs. 11 and 12. For p< 1-4 or 1°5 there is quite good agreement, 
but further out Michelson’s intensities are higher than ours and his intensity 
gradient appreciably less. 

A relative photometry of the corona at the same eclipse was also carried out 
by G. M. Nikolsky (44), again at a site near Cheeli. He worked with a focal 
length 195 cm and in the wave-length region A 4600-6600, which is not quite the 
same as ours. Unfortunately, the original plate calibration, with a tube sensito- 
meter, was not successful and a substitute calibration had to be taken from 
another instrument, but with the same plate and processing. Nikolsky’s paper 
gives two sets of isophotes from two plates, covering 1°25 <p< 312, and after 
adjusting his scale to ours the results are summarized in our Fig. 21. His gradient 
is less steep than ours, the difference at both ends of our common range amounting 
to about +0-2 in lg(K+F). 
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Summary 


Signal-to-noise equivalent numbers of stored photons and equivalent 
quantum-efficiencies are derived for four emulsion batches from homogeneous 
sensitometry and granularity data. ‘The highest monochromatic efficiencies 
found are o’g per cent in the blue and 0-6 per cent in the red. Peak efficiencies 
occur at densities of only o-1 to 0:2 above fog. The peak equivalent number 
of photons stored ranges from 0-03 per jp” for the “ fastest’? emulsion in the 
blue to 0-4 per p? for the finest grained material in the red. The distinction 
between exposure for maximum efficiency and for maximum equivalent 
storage is emphasised. 





1. Introduction.—Considerable interest has recently been shown (1) in 
comparisons between different types of radiation detector, and in particular the 
importance has been recognized of comparing the sensitivities of the newer 
kinds of photoelectric image-detectors with those of photographic emulsions. 
In the early development of a new device it is common to make empirical 
comparisons with the older devices that it is intended to supplement or replace ; 
for example, in the early development of mechanical power, pulling matches 
were arranged between traction engines and horses, and similarly between 
paddle- and screw-driven steamers. ‘The limitations of these trials are obvious 
to us now. If conducted by ‘‘tug-of-war’’ they really measured force and not 
power, and they were too specific to the particular circumstances of the trial to 
allow account to be taken of the variation of performance with the conditions 
of working. For example, one engine might develop more power than another 
when both were running fast on a light load, while the second might be the more 
powerful when they were running slowly on a heavy load. A valid and compre- 
hensive comparison was possible only when an adequate definition of what is 
meant by ‘‘ power” had been formulated, and when a theory was available 
which indicated what quantities must be measured, and how the measured values 
must be combined, in order to determine the power in the defined sense as a 
function of the running conditions. 

Empirical comparisons between radiation detectors are subject to similar 
limitations, and in consequence it may happen (as anyone knows who has carried 
out such trials) that a marked apparent superiority of one detector over another 
in one test may be reversed in another test under different conditions. Before 
adequate comparisons can be made, it is necessary first to define what is meant 
by ‘‘sensitivity’’. Procedures then follow for abstracting ‘‘sensitivity’’ in 


the defined sense from quantities which can be more directly measured, just as 
in the case of mechanical prime-movers it could be shown that the product of 
force and velocity was the appropriate measure of ‘‘ power’’. 
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In any new technology, measurements may for a time be made in arbitrary 
units (such as the jar for electrical capacity, the inch of mercury for atmospheric 
pressure or the candle for light) before the link is made to corresponding absolute 
physical units (farad, bar and photon sec™! respectively). Watt’s (2) calibration 
and definition of the horse-power in foot-pound units represents this step in the 
measurement of power, and it is commemorated by the use of the name watt 
for the practical c.g.s. unit of power. ‘The assessment of radiation detectors has 
followed this pattern in so far as the quantities commonly used as measures of 
sensitivity have been mainly those which give a rough idea of the engineering 
performance under restricted conditions of use. For example, lumen respon- 
sivities do this with respect to use with tungsten lamp illumination. Such 
measures take no direct account of noise, and accordingly they are not sufficiently 
realistic for present-day applications where (with few exceptions) the amplification 
that can be used with a detector is limited only by the noise level at its output. 
Thus amps/lumen in a CdS photoconductive cell can represent a worse faint-limit 
than »A/lumen in a photoemissive cell. A further disadvantage of these measures 
is illustrated by the fact that Sb-Cs and Ag—Cs—O photocells have similar 
pA/lumen ratings, but their performances on an average star are strikingly 
different. ‘The reason is of course that lumen responsivities are inapplicable 
to any other illuminant than tungsten lamps. In the absence of more realistic 
and general measures of sensitivity, direct empirical comparisons may be the 
only expedient available for relative assessment of detectors. But as soon as it 
is feasible to make the necessary measurements in absolute physical units, curves 
of sensitivity against the relevant parameters of use become the most informative 
means of comparing different detectors. Direct comparisons, if carefully 
controlled and interpreted, can at this stage serve as valuable spot-checks on the 
accuracy of the sensitivity assessments. 

It is not obvious at first sight how the sensitivity of a radiation detector 
can be most suitably defined. Evidently no detector, however nearly perfect, 
can give more than a certain signal-to-noise ratio in an image, or in an 
intensity-measurement, at a given light level, because its performance is limited 
by the spontaneous fluctuation in the relevant photon-rates. A quantity, first 
formulated in the present connection by Rose (3), takes specific account of this 
circumstance and appears to have outstanding advantages as a measure of 
sensitivity. This quantity, which I have termed ‘‘ equivalent quantum-efficiency ”’ 
in previous discussions of some of its astronomical applications (4, 5), may be 
thought of as the ratio of the light level at which an ideal detector would yield a 
given quality of image or accuracy of measurement, to that which the detector 
under consideration actually needs in order to give the same performance. ‘This 
ratio lies in the range 10~!-10~% for most modern detectors of radiation. A more 
precise definition is given in Section 2, followed by discussion of its application 
to photographic emulsions. 

It turns out that the equivalent quantum-efficiency of an emulsion at a given 
photographic density can be calculated from the exposure in ergs/cm? or 
photons/cm? necessary to produce this density, the granularity at this density, 
and the slope (gamma) of the ordinary characteristic curve at this point. Of these 
quantities, the granularity is the least accessible. Its measurement requires 
rather specialized equipment, and, as Zweig (6) has remarked, a rather large 
number of measurements across the plate must be made in order to attain 
sufficient statistical accuracy. These and other circumstances have deterred 
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observatories from attempting granularity measurements, and as the manufacturers 
do not quote granularity quantitatively in their published characteristics of 
emulsions, the data necessary for a straightforward determination of photographic 
equivalent quantum-efficiencies have not hitherto been available. Consequently, 
the estimates that have been made have been indirect, or have involved combining 
sensitometry and granularity data that were not homogeneous (in the sense 
that they did not refer to identical emulsion batches and processing), or have 
depended on incomplete data, so that assumptions that may prove to be inaccurate 
have had to be introduced into the interpretation. In a previous paper (7), 
I have made estimates in these ways for a number of commercial and astro-type 
emulsions, and presented the results as curves of equivalent quantum-efficiency 
against exposure time in order to show the effect of reciprocity failure. These 
curves indicate that the peak performances of the commercial and astro-type 
emulsions for which data were then available are not very different, though the 
peak of course occurs at very different exposure times for the two types. The 
relative values of equivalent quantum-efficiency shown in these curves are 
believed to be correct within a factor two, but the absolute values may be less 
reliable. ‘This is the more unfortunate since, as I pointed out earlier (4), there 
has been some conflict of evidence as to whether the best equivalent quantum- 
efficiencies for photographic emulsions are around o-r per cent or I per cent. 
Recently, the consensus of opinion seems to have favoured the lower value (1). 
In my own work, however, I have found increasing evidence in favour of the higher 
one, and have presented some possible reasons for the discrepancy (7). One 
such reason is that peak monochromatic values are needed for comparison with 
the photoelectric values usually quoted, whereas some workers have used 
‘‘ white-light ’’ values. 

However, the evidence I have previously adduced in favour of a value of about 
1 per cent for the equivalent quantum-efficiency has been subject to the 
difficulties of indirect, inhomogeneous or incomplete data referred to above. 
Hence, the individual estimates should certainly be treated with reserve: and, 
although in aggregate the persistence with which they yielded values of 1 per cent 
or more may seem significant, nevertheless there has remained ground for 
reasonable doubt as to the absolute performance of photographic emulsions. 

The Eastman Kodak Company has now kindly made available homogeneous 
granularity and absolute sensitometry data on named emulsion batches. ‘These 
new data allow performance estimates for emulsions to be made on a much more 
satisfactory basis than hitherto, and the estimates obtained from them form the 
substance of the present paper. It is natural that the present data, the first of 
their kind to be published, should refer to commercial emulsions, which are the 
makers’ principal concern. Although much astronomical photography (for 
example, photography of planets and of solar phenomena) uses these emulsions, 
and their use may increase if cooling is adopted in order to control reciprocity 
failure, nevertheless the majority of astronomical photography involves long 
exposures and the use of emulsions with low reciprocity failure suited to these 
exposures. Itis very desirable that the present good beginning in the measurement 
of the quantities needed to determine equivalent quantum-efficiency should be 
extended to astro-emulsions of this type. The potential saving in observing time, 
resulting from a detailed evaluation of the optimum conditions of exposure 
(see Section 4), would fully justify such a programme. 
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rers 2. Definitions.—If N photons are in principle necessary to obtain a given 
: of measurement, and if an actual detector in fact requires (with the same informa- 
hic tional coding) O incident photons, then the equivalent quantum-efficiency €¢ is 
tly, defined as 
‘ing €e=N/Q. (1.1) 
nse More specifically, suppose that QO photons are available to a detector and that 
a they yield a single intensity measurement with ratio R of signal to r.m.s. noise. 
re The number of photons necessary in principle to yield this accuracy N= R?, 
(7), and therefore, by the definition (1.1), the equivalent quantum-efficiency is 
e 

a <o= R2/0. (1.2) 
ese The meaning of this relation may be illustrated by comparing the performance 
ype of the actual detector with that of a hypothetical detector which acts on the O 
the photons in such a way that a fraction, ¢«, of them are fully used and the rest wholly 
he lost, so that the number of photons effectively used is eO and the signal-to-noise 
are ratio R, is given by R,?=«Q. If the fraction is chosen equal to the equivalent 
less quantum-efficiency €., we obtain R,?=«Q=e«,.Q = R’, and hence the performance 
ere of the hypothetical detector is the same as that actually attained. The quantity 
m- N=eeO may be called the equivalent number of photons used or ‘“‘stored’’ in 
nt. the measurement; it provided a convenient measure of the quality of the 
I). measurement. 

her Equivalent quantum-efficiency is in general a function of conditions of use. 
ne For example, its variation with total exposure, exposure time, wave-length of 
ith light, and spatial wave-number of the distribution sought are concerned 
sed respectively with the overload characteristics, reciprocity failure, spectral 

sensitivity and spatial resolution of the detector. 

out Equivalent quantum-efficiency «g and equivalent number N of photons 
the contributing to a measurement are, as mentioned in Section 1, concerned with 
ve. a rather different aspect of photography from ordinary sensitometry, and one that 
id, is of more practical importance for astronomical measurements. Sensitometry 
ont is concerned with the average relationship between light intensity and the resultant 
for photographic density. The quantities «, and N are concerned with how good a 
ns. measurement (in this case in the r.m.s. sense) can be obtained, in the presence 
Us of the fluctuations due to granularity, when systematic distortions introduced 
se by the non-linear relationship of density to exposure have been removed, as they 
re can and should be in astronomical measurement by means of calibration exposures. 
he It is to be expected that both «,g and N will have maxima for some particular 
of exposure; in ordinary photographic parlance ‘‘under-exposure’’ corresponds 
he roughly to exposure less than that which gives maximum N, and “‘ over-exposure ”’ 
‘or to greater exposures than this optimum. Evidently it is a non-sequitur from the 
1S, definition (1.1) to identify N with the number of grains made developable. 
ity (On certain assumptions (7) N is approximately equal to y?/D? times the number 
ng of developable grains. ) 

Se The formulation which follows in Section 3 is expressed in terms of approximate 
nt differential equations and is applicable as it stands to low-contrast images. Since 
be both e.g and N depend on exposure, the corresponding finite-difference equations 
1e, should be used in the discussion of high-contrast images. For any finite contrast 
re there is an exposure below which e, and N are effectively zero because no 


appreciable photographic effect is produced, and an exposure above which they 
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are effectively zero because the emulsion is everywhere fogged to maximum 
density. 

3. The data and calculations.—The data refer to four emulsion types which 
are identified in the notes to Table I. The sensitometric data consist of curves 
of density D9 in log, units against log,) exposure FE in ergs/cm?. These curves 
are given for quasi-monochromatic exposures at two effective wave-lengths, 
A=430 and 607mp. They refer to 15 sec exposure, but factors for correcting 
to 1/10 sec (approximately the optimum exposure) are given on the curves, and 
have been applied in working out the ¢g. Thus the e, values refer to nearly 
optimum exposure times. 

The granularity data give the r.m.s. fluctuation in density AD,, for various 
densities D,, and for circular scanning apertures of diameter a varying from 10 
to4op. Itisstated that AD,, isa function of D,, that does not depend appreciably 
on the wave-length of the exposing light. 

Let the slope of the sensitometric curve be written 


_ GdDy — dD _ EdD 

Y~ dlog,£ dinE dE ea 
where D is the density measured in log, units. Then if the r.m.s. density 
fluctuation is AD for a scanning area of A square microns, the signal-to-noise 


ratio R of the measurement of intensity over area A is given (approximately) 


by 





EE dD_y 


= XE ~ XD‘ dE ~ AD is 
and by equation (1.2), the equivalent number N of photons stored per square 
micron is 
_ tfy\ 
N=-=(—}. 2. 
x(z5) (2:3) 


In terms of the density fluctuation expressed in log,) units, and of the diameter a 
of the scanning circle expressed in microns, this formula takes the form 
r_ {o49y)? 
N= ¢—*7"\. . 

la AD jp (2-4) 
The equivalent quantum-efficiency is found by dividing N by the number of 
incident photons per square micron. Specifically, with the same units as in 
equation (2.4), exposure E expressed in ergs/cm?, and A in microns 


' = 2 7} 
€e=4°77 X 10 (@AD,,)°EA | 
2 
= 111x104 at A=0'430pn, f (2.5) 
10 
= 9°Q X 10 a iE at A=o-607p. | 
10 } 


The results are presented in Table I and Fig. 1. In both presentations, the 
emulsion type is identified in a mnemonic manner, and the exact identifications 
are given in notes to the table. The characteristics of these emulsions should 
not be taken as applicable to other batches of nominally the same type. The 
manufacturers emphasize that while the data are representative of the emulsions 
at the time of manufacture it must be recognized that the characteristics of 
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Type A(p) 
Royal-X 0°430 
0607 

Tri-X 0°430 
0°607 

Plus—X 0°430 
0°607 

Pan-X 0°430 
0607 


(3) 
a(p) 
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O|p? 
4°5 
24°2 
119 
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39°3 
164 
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134 


807 


(11) 
51 


165 


718 


(12°2) 


272 
1,460 


(31) 

119 

473 
2,270 


(26°6) 
116 
533 


25440 
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D(log 10) 


1°36 


0°28 
0°58 
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0°28 
0°52 
1°02 
1°62 
028 
0°52 
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1°62 
0°26 
0°48 
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(7) 


e,(per cent) 


o’71 
0'073 
o'o12 


(0°12) 
0°59 
0°085 
o’ol1o 


(0°66) 
O31 
O21 
O14 
0°057 
0°037 
0°034 
0°0052 
0°0052 
0°0025 


(0°02) 
0°40 
0°27 
O'l7 
o’ogI 
0°059 
0°055 
0:0068 
0°0068 
0°0035 


(0°22) 
o-41 
lobe fe) 
0°024 


(0°25) 
0°29 
0076 
O'O15 


(0-39) 

(0°28) 
0°20 
0°16 
0°060 
o'051 
00065 
0:0086 


(0094) 

(0069) 
O'13 
o'10 
0°047 
0°040 
0°0055 
0°0075 
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lotes to Table I 


Col. 


1 Emulsion Type. Royal-X=Kodak Royal—X Pan Film, Code 6128. 
Tri-X=Eastman Tri-X Panchromatic Negative Film, 
Type 5233. 
Plus-X=Eastman Plus—X Panchromatic Negative 
Film, Type 4231. 
Pan-X = Kodak Panatomic—X Film, Code 5240. 
All films manufactured February 1957. 


2 Wavelength A of emposing light in microns. 

3. Diameter a of scanning aperture in microns. 

4 Exposure Q in photons per square micron. 

5 Photographic density Dy, in log, units. 

6 Equivalent number N of photons stored per square micron. 
7 Equivalent quantum-efficiency ¢, as a percentage. 


products of the same name may vary within manufacturing tolerances and may 
change significantly as improvements are incorporated in them. 

In Table I, the equivalent quantum-efficiency €, (per cent) and the equivalent 
number of stored photons N (per sq. micron) are listed for the densities at which 
the granularity AD,, has been measured. It seems possible that the smallest 
of these densities is, for each emulsion, intended to refer to the fog level. However, 
most of the D,,—log E curves appear to show a small slope at these densities, 
and the corresponding N and €, values have been included in parentheses in the 
table. The values of N and eg would be independent of the size of the scanning 
aperture for an emulsion obeying the Selwyn law of granularity ; for real emulsions 
there may be a size of aperture which maximizes Nande,. The practical usefulness 
of this optimum depends on the extent to which image detail is blurred by 
turbidity. An adequate discussion of this effect requires a knowledge of the 
spatial spectra of the photographic response and of the granularity fluctuations, 
and cannot be attempted with the data at present available. ‘The values of N 
and €, for aperture diameters a= 10, 20, and 40 are given separately in Table I 
whenever, on internal evidence, there seemed to be a significant dependence 
ona; in the remaining cases the mean values over all available a has been given. 
In particular the values for the Royal-X type emulsion show no systematic trend 
with aperture size. 

In plotting the curves shown in Fig. 1, it has been necessary to make an 
interpolation between the densities for which AD, values are available. It is 
clearly not satisfactory to interpolate the N and e, values directly, owing to their 
relatively rapid change with exposure. It was found, however, that the values 
of AD,,/a could be interpolated without significant uncertainty. The available 
Dy. and y values were then used to calculate N and €, at intervals of 0-2 log) units 
in E, All the curves refer to the mean values over all the available aperture sizes. 
This reduces the effect on the curves of the statistical fluctuation in the AD, 
values; in particular, it helps to prevent these fluctuations from causing 
spuriously high peak e, values to appear. 

4. Discussion.—The results support the view that the equivalent quantum 
efficiencies of the best emulsions can be near to 1 per cent. If this conclusion 
is true, it implies that it will be much more difficult for other kinds of image- 
detector to yield sensitivities superior to that of photography than it would have 
been if the lower value more commonly assumed, 0-1 per cent=€e, had been 
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ig Fic. 1.—Photographic density D in log, 9 units, equivalent number N of stored photons per 
square micron, and equivalent quantum efficiency €, as functions of log, exposure. The four 
rows refer to four emulsion types which are identified in the notes to Table I. The left-hand set 

m | of curves refer to exposure to light of wave-length 430mp, and the right-hand set to607mp. Note 

mn that the scale of the D and ¢, ordinates is the same in each panel, but that of Nis not. The abscissa 

a represents log, exposure in ergs|cm* to the same scale as the ordinate D. The zero-point of the 
abscissa varies from panel to panel; the curves are most conveniently interpreted by reference to D 

fe as parameter. 
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correct. ‘This in turn means that the necessary development of new detectors, 
such as photoelectric image tubes, will probably take longer. In addition, the 
gain in sensitivity over photography that can be expected from presently fore- 
seeable technical developments is by a factor of tens rather than hundreds. This 
gain in sensitivity, combined with the other advantages of photoelectric methods, 
still appears sufficient to make the development of photoelectric image-detectors 
for astronomy one of the most important advances in observational methods 
being undertaken at the present time. 

The curves confirm previous suggestions (7) that the dependence of ¢€. on 
exposure may be very sharp, and that the maximum €e occurs for rather low 
photographic densities, usually only o-1-o-2log,, units above fog. This 
emphasizes the advantage, where possible, of coding the information presented 
to the emulsion in such a way as to reduce the contrast range in brightness, 
For example, an extended sinusoidal pattern can be more efficiently photographed 
than a ‘‘point’’ image carrying the same information. 

In so far as the maximum of N corresponds to the negative-exposure which 
results in the least grainy print of given contrast, these conclusions accord well 
with ordinary photographic experience. As Mees (8, pp. 866-8, esp. Section 4) 
has pointed out, the optimum exposure in this sense corresponds to a thin negative 
and he accordingly recommends the minimum negative exposure that is 
satisfactory on other grounds. Print quality as ordinarily judged involves 
grey-scale reproduction as well as graininess. ‘The minimum exposure which will 
give a subjectively satisfactory print depends on the contrast range of the image, 
one condition being that the darker areas should be reproduced with adequate 
signal-to-noise ratio. There is a range of exposures greater than this minimum 
over which print quality is consistently judged as ‘‘excellent’’ (Mees (8), 
pp. 766-721), and in ordinary photographic practice, where length of exposure 
is not expensive, the aim is rightly to give an exposure corresponding to the 
middle of this plateau. ‘This ensures that moderate errors in exposure do not 
degrade the print quality appreciably. The efficiency in this region, however, 
falls at least as fast as the reciprocal of the exposure, and this loss in efficiency 
should not be tolerated in astronomical photography. Some modification of the 
usual techniques of measurement by microphotometer may be required in order 
to exploit the higher efficiency occurring at low densities. 

The curves also confirm (7) that there is often, especially for the finer-grained 
materials, a marked difference between the conditions leading to maximum 
storage and maximum efficiency. At the exposure which gives greatest storage, 
the efficiency is commonly 1/5 to 1/10 of its maximum value. Thus exposing 
a plate until with further exposure no new star images appear does not make the 
most efficient use of the available observing time, and €, values based on such 
exposures may be misleadingly low. 

During the final revision of this paper, I received MSS by Dr R. Clark Jones 
(dated 1958 January 1; 1957 December 1; and 1958 March 4) in which he 
analyses and compares the performance of photographic negatives, the human 
eye, and image-orthicon camera tubes, in terms of what he calls ‘‘ detective 
quantum-efficiency ’’, which is identical with the present ‘‘ equivalent quantum- 
efficiency’’. He concludes that human vision has an equivalent efficiency of 
I per cent at 1 ft lambert, and that the camera tube may attain an efficiency of 
2-4 per cent under optimum conditions. The rough estimates which I previously 
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put forward (7) accord fairly well with these values. Jones’s analysis of photo- 
graphic performance is based on the same data as are used here (but for the blue 
only), and he derives a working formula differing in form but identical in content 
with the present equations (2.5). Apart from some small differences, due to 
his use of a smoothing process which I omit, he derives the same values of 
equivalent efficiency as are presented here, and draws similar conclusions from 
them. He regards the background density D as being due to a pre-exposure 
upon which the signal giving the change in density AD is superposed, whereas 
I regard both D and AD as being due to a single exposure to a low-contrast signal. 

I take pleasure in acknowledging my indebtedness to the Eastman Kodak 
Company for making available the data on which this investigation is based, and 
to Dr J. A. Leermakers and Dr O. Sandvic for their help and advice. It is also 
a pleasure to acknowledge my debt to Dr R. Clark Jones for the privilege of 
seeing his MSS in advance of publication, and for much helpful discussion. 


The Observatories, 
Madingley Road, 
Cambridge: 
1958 March 12. 
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TRANSITION PROBABILITIES FOR FORBIDDEN LINES 
OF Ni II AND Ni III 


R. H. Garstang 


(Communicated by the Director of the University of London Observatory) 
(Received 1958 April 8) 


Summary 


Transition probabilities have been computed for magnetic dipole and 
electric quadrupole radiation for all transitions between levels of the 3d° 
configuration of Ni 111 and between the levels of the 3d*4s and 3d°® configura- 
tions of Ni II. 





1. Introduction.—In a previous paper (4) the first results were given of a 
programme of calculations of transition probabilities for forbidden lines of 
astrophysical interest arising from atoms in the first long period, data being 
obtained for [Fe 111] and [Fe v]. Pasternack (13) had earlier presented data 
for [Cr 1v], [Mn v], [Fe vi] and [Fe vir]. In the present paper results are given 
for transitions in [Ni 1] and [Nii]. The forbidden lines in Ni 11 occur within 
and between the configurations 3d*4s and 3d°, those in Ni 111 within the con- 
figuration 3d’. The calculations on Ni 111 were originally undertaken to assist 
in the identification of lines of this ion (following upon the search by Swensson 
(2x)). Ni 1 was studied as a simple example of an ion for which transitions 
of the form d”—d"~'s and d"-1s—d"~—s have to be considered, and for which 
configuration interaction should be included in the calculations. Forbidden 
lines of Ni 11 occur in a number of celestial spectra. 

2. Calculations on Ni III.—The transitions occur between the levels of the 
3d® configuration. The 3d’4s configuration is well separated from the 3d* 
configuration, and the former configuration was neglected in the calculations. 
The spin-orbit matrix is the same as that for the 3d? configuration given by 
Condon and Shortley (2, p. 269) except for the reversal of the sign of £,, (2, p. 299). 
The electrostatic energy matrix is diagonal, and the diagonal elements E(°F), 
E(?}D), EP) and E(!G) were treated as independent parameters, to be deter- 
mined along with ¢, to give the best fit of calculated and observed spectroscopic 
energy levels. ‘Theoretical formulae for the electrostatic energies (2, pp. 202 and 
299) were used to obtain an estimate of E('S). The observed energy levels 
were taken from Shehstone (18). The values of the parameters finally adopted 
are listed in Table I and the resulting energy levels are given in Table II. After 
completing the fitting of the observed energy levels to the theory, the transforma- 
tion matrix which diagonalizes the complete energy matrix was determined. 
The matrix of the square roots of magnetic dipole line strengths in LS-coupling 
was set up using formulae given by Shortley (20). The corresponding matrix 
for electric quadrupole radiation was given by Pasternack (13, with two corrections 
noted in (3) and a change in all the phases for the almost-closed shell). The 
transformations to intermediate coupling were performed and the resulting 
line strengths converted into transition probabilities in the standard manner. 
In this last stage of the computations a value of the radial quadrupole integral 


+= | % r2P2(3d)dr 
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Adopted parameters for Ni 111 (3d®)(in cm-) 
E@F) 1022 
E(:D) 14365 
E(P) 16610 
E('G) 23089 
E('S) (64500) 
Sa 75 
Tasce II 
Observed and calculated energy levels for Ni 111 
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was required. This was estimated by calculating s, from wave functions given 
by Hartree for V*#(8), Mn*#(7) and Zn+*(8), and interpolating (s,"/r), where 
r is the mean radius of the 3d wave function as tabulated by Hartree (8, Table I1), 
The estimate obtained in this way for Nii iss,=1-10. The final values of the 
magnetic dipole and electric quadrupole transition probabilities are given in 
Table III; they are spontaneous emission probabilities in reciprocal seconds. 

3. Calculations on Ni [J.—The transitions of interest are between and within 
the configurations 3d° and 3d*4s. ‘The method used was essentially similar 
to that for Nii. ‘The spin-orbit matrix for the d*s configuration was obtained 
from a formula given by Trees (25, footnote g). The matrix for d*s had been 
given by Marvin (10), but his phases, though mutually consistent, are not in 
accordance with the standard choice of Racah, and the importance of correct 
phases in line strength calculations necessitated the recomputation. It was 
decided to allow for configuration interaction between d® a*D and d*s b?D. The 
matrix element of this interaction was obtained by Ufford (26); it was verified 
as to both magnitude and phase by the writer. There is no spin-orbit interaction 
between the d® and d*s configurations. ‘The electrostatic energy matrix is 
diagonal except for the configuration interaction element, and the diagonal 
elements were treated as independent parameters to be determined along with 
C, to give the best fit of calculated and observed spectroscopic energy levels. 
The observed energies are taken from (12). The highest level, d*s 7S, is not 
observed, and its position was predicted using theoretical formulae (15, and 2, 
p- 203). 

The determination of the configuration interaction parameter (H, in the 
notation of Racah (14)) was rather difficult. However, values of H, were found 
in the literature for Ti 11 (9), V 11 (g), Cr 1 (17, 25) and Mn (25). It was 
noticed that H,F,* was approximately constant along this sequence of ions, F, 
being the usual Slater parameter; this was used to predict H, for Nit. The 
value adopted was H,=—55cm~'. The sign is important: the phase of H, 
must be consistent with the phases of the transition integrals. ‘The remaining 
parameters were adopted by the usual process of trial and error; the final values 
of the parameters are given in Table IV. The observed and calculated energy 


TaBLeE IV 
Adopted parameters in Ni 11 
E(?G) 32492 
E(‘F) 9418 
E(?F) 14142 
E(b*?D) 24037 
E(‘P) 24368 
E(?P) 29154 
E(?S) 63900 
E(a*D) 612 
C (d's) 661 
C ,(d°) 603 
H, =§5 


levels are compared in Table V. The fit is quite good. As a further check 
on the calculations, Landé g-factors were calculated, and are also given in Table V. 
Although not perfect, the calculated g-factors are much nearer to the observed 
values for the critical b?D and ‘P levels than to the LS-coupling values. (‘These 
levels have been discussed by Racah (15) and by Shenstone and Wilets (19); 
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they suggest that the assignments of the 6D,, and ‘P,, levels should be inter- 
changed. Our calculations confirm this if one assigns to a mixed level the 
designation of the largest constituent LS-coupling level. However, in order 
to avoid confusion, the level designations in this paper have been kept as in 
Atomic Energy Levels (12) and as in the Revised Multiplet Table (11)). 
TABLE V 
Observed and calculated energy levels and Landé g-factors 


Landé g-factors 


= Energy Levels tail : a 
Term J Cilicia Cileaieied Observed LS- Calculated 
Coupling 
aD 23 ee ° —— sist 1'200 1200 
1} 1507 staid 1508 — a 0800 ~— 0800 
‘F 4} 8304 _ 4, ©6840 — 1°35 1333-1333 
3} 9331 __ hed 9332 _ 81 1°24 1'238 1°236 
2} 10116 rte | 10113 _ 7 3 1023 1'029 1°027 
14 10664 54! 10661 54! 0°397 0°400 07402 
I 34 13550 13546 1°14! 1°143 1°144 
KB Je “a ¢ Je ~_ < 
24 14996 — 14997 1451 0866 0°857 0861 
#D 2h 23108 _ 6gg 23139 ~~ 1°428 1*200 1°409 
1} 23796 m 23752 3 1'045 0°800 1005 
‘p j 24836 P 24873 ee 2°667 2667 2°664 
14 24788 pos 24800 oa 1°498 1°733 1°540 
2} 25036 25047 1°368 1'600 1°388 
i og 1} 29071 29094 1°322 1°333 1°318 
} 29593 ae 290574 — 0670 ©=—.0'667 ~— 0669 
*G 4h 32499 4, 32510 — 1135 II! Ir! 
33 32523 32513 0895 0889 0889 
2S h ea (63967) Sef 2000 2000 


The transformation matrix which diagonalizes the energy matrix was 
determined numerically. The matrices of the square roots of line strengths 
were obtained for magnetic dipole (20) and electric quadrupole (3, 5) radiation ; 
they were transformed to intermediate coupling and transition probabilities 
computed in the standard manner. The chief difficulty was the estimation 
of the quadrupole transition integrals. So little information on the behaviour 
of the wave-functions of 4s electrons was available that the transition integral 
for 3d — 4s transitions could not be estimated with adequate accuracy on the basis 
of results for other atoms. Fortunately, D. F. Mayers has recently completed 
programming self-consistent field (with exchange) calculations for an electronic 
computer, and he very generously offered to compute wave functions for the 
34° and 3d*4s configurations of Nit. Using his wave functions it is found that 


t= f ” »2P2(4d)dr = 1-42, for 3d in 34°, 
0 

5,= . r?P?(3d)dr = 1-19, for 3d in 3d*4s, 
0 


Sg= {- r?P(3d)P(4s)dr = — 2:26, for 3d in 3d° and 4s in 3d%4s. 
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These values were adopted in our calculations. The final magnetic dipole 
and electric quadrupole transition probabilities are given in Table VI; the 
results are spontaneous emission probabilities in reciprocal seconds. Under 
the multiplet designation we have inserted the number in the Revised Multiplet 
Table (11). 

It seemed desirable to check the configuration interaction parameter H, by 
calculating its value using Mayers’ wave functions. The result obtained was 
H,=-—70cm~1. The agreement with the empirical value is as good as can be 
expected. The effect of including H, was examined during the calculations 
It was found that its inclusion had little effect on transitions within the d*s con- 
figuration, but transitions between the d*s and d® configurations were affected 
to a substantial extent; neglect of H, would have given rather unsatisfactory 
results for the latter transitions. Configuration interaction had only a slight 
effect on the calculated energy levels. 

4. Comparison with observations.—Swensson (21) made a search for lines 
of [Ni 111] in various celestial objects. He found some evidence for the lines 
3F,-8P, and °F,-*P,, but the results were inconclusive. Osterbrock (private 
communication) has renewed the search for these lines, and finds a line at 
A6401°6 in the Orion nebula which may be *F,-*P,; he found no evidence of 
3F,-8P,. The transition probabilities in Table III lead us to think that these 
two transitions might appear weakly in objects with suitable excitation, and the 
identification of *F,—°P, is a plausible one. The lines 1G,-°F, (A4326-2) and 
1G,-5F, (A4596°8) in the blue region of the spectrum should be the most easily 


TABLE VI 
Transition probabilities of [Ni 11] 
Transition A,, Ie Transition A. Ae 
a’D-a®D 23-14 = o055 1°4X 1077 *F-"F) —s_ 34-24 = 01047 6°3 x 107° 
2T>-4 ae cee -6 
“Stine ee ee lee 
2732 2 (3F) 14-2} 9:5x10-? 0-090 
13-33 se 9°9 x 107° a ee ' 
24-24 =393-7x 10° )=— 14 x 1078 — . oe _—— 
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24-15 =rrgx 1075 =2-2 x 107% ‘F-hD 44-2} Avs o-019 
th-13 gx 1075 = 3-3. x 1075 (7F) 34-23 0°17 4°2 x 107% 
‘F-“F = 444-34 = 0027 1°5 Xx 107° 29-24 OOI2 8-4 x 107 
443-24 a 3°0 x 107° 13-25 95x 10-% gg: x 10 
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34-1} ae 1-7 X 1079 24-14 = 0°39 2°3 x 107 
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aD-*F 24-34 2:9x1074) 014 °F-b7D 34-24 ~=0'078 zaxro~° 
(2F) 14-34 ii 7°6 x 1073 (11F) 24-24 0084 4'2 x 1078 
24-24 26x1074 0:062 34-14 ea 9°3 x 107° 
14-24 3=5:5x 1078 ovr! 24-14 = 0058 3°8 x 1075 
4F#2F 43-31 0086 ; I'4X 10"? bD-b2D 24-1} 53x 10-9 ~— 1-3. x 107 
34-34 5:0x107% 2:0x 107 
23-34 0-014 3°9 x 107° a@D-'P 2}-24 0-42 0°42 
14-34 9°5 x 10719 (4F) 14-24 76x107% 014 
4}-24 ae 1°8 x 1078 24-14 0034 0°042 
34-24 3:4x107 2:1 x 1078 rh-14 ov o'14 
24-24 3 37:5x1073 4:3x1078 2h- 4 oe 3°9x 107° 
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observable lines of [Ni 1]. Their transition probabilities (Table III) are high. 
There is a line in 7 Carinae and RR Telescopii (24) near 44596; this can now be 
identified with 1G,-°F;._ The line 1G,—°F, is almost certainly blended with a 
- strong line of [Ni 11] and the [Ni 111] line is unlikely to be observable in other 
j 
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objects with strong [Nit] lines. The lines 1D,-°F, (A8500) and 4D,-*F 3 (A7890) 
should be strong, but they lie in an unfavourable region of the spectrum and do 
not appear to have been observed. 

Swings (22) first identified [Ni 11] lines in stellar spectra. Some lines in the 
multiplets a*D-a?P, a*D-a?G and a?D-6?D were found. He also predicted 
that the infra-red multiplet a?D-a?F should be strong. The latter multiplet 
was observed in great strength in » Carinae by Thackeray (24). A number of 
other observers have identified one or more multiplets in various peculiar stars, 
including » Car (6, 24), Boss 1985 (1, 22, 23) and WY Vel (16). There is fair 
agreement between the relative strengths of lines of various multiplets observed 
in these stars and the corresponding theoretical strengths. There are, however, 
some discrepancies: the two most notable are the lines a?D,,—°F,, and 
a’D,.,-b*D,,, for which the observed lines are stronger than would be expected 
on the basis of the theoretical transition probabilities. It seems probable that 
the theoretical values are too small. 

Acknowledgments.—The writer is greatly indebted to Mr D. F. Mayers for 
undertaking the calculation of wave functions for ionized nickel and for com- 
municating his results in advance of publication. 
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Summary 


A discussion is given of methods which may be used to compute Gaunt 
factors for free-free transitions in the neighbourhood of any positive ion. 
The equivalence of the formulae of Sommerfeld and of Menzel and Pekeris 
for hydrogen-like systems is demonstrated, and a number of approximations 
to them discussed. It is shown that the rapid computation of numerical 
values of monochromatic Gaunt factors for any electron velocity is assured 
by means of an extensive table of values supplemented by simple asymptotic 
expansions. 

The effect of screening on Gaunt factors is estimated for long wave-lengths 
by a simple formula derived in the Born approximation with the aid of an 
atomic potential function of the form r~!ZXa, exp (—f,7r). The method is 
illustrated by examples in which Thomas—Fermi and self-consistent potential 
functions are expressed in this form. 

A simple procedure for calculating Gaunt factors averaged over a 
Maxwellian velocity distribution is given using a quadrature formula based 
on evaluation of the integrand at the zeros of Laguerre polynomials. Three 
integration points are found to be adequate to give accuracy comparable to 
that attained by more elaborate methods. 





Introduction 

The purpose of this paper is fourfold. In the first place, we wish to examine 
the formulae available for computing the continuous absorption of radiation by 
free-free electron transitions in the neighbourhood of an atomic nucleus. 
Secondly, we shall uncover the connection between the formulae of Sommerfeld 
(6) and Menzel and Pekeris (5). This connection has been obscure for many 
years, and, in view of the competing claims of the two expressions in calculations 
of free-free absorption and emission, its elucidation seems overdue. Thirdly, 
the role of simple approximate formulae is discussed, paying particular attention 
to the range of photon frequency and electron velocity in which they are valid. 
The effect on the free-free cross-section cf the screening of the nucleus by the 
surrounding electrons is of particular importance for long wave-lengths, and 
approximate formulae are developed by which order of magnitude estimates 
can be made in this limit. Finally, it is hoped that the tables presented here 
together with the simple approximations discussed will be sufficient to enable 
anyone wishing to compute a cross-section to do so rapidly and with ease. 

In addition to the problem of obtaining cross-sections for specific frequencies 
and electron velocities, many astrophysical problems require an average to be 
taken over a statistical distribution of electron velocities. A simple method of 
doing this for the Maxwell distribution has been developed, and the final section 
of this paper is devoted to explaining this technique and to a few numerical 
examples chosen to illustrate the methods described earlier. 








242 I. P. Grant, Calculation of Gaunt factors Vol. 118 


The first calculations of the cross-section for free-free emission were made 
by Kramers (1) and Wentzel (2) in terms of a classical description of the particle 
orbit. The cross-section for spontaneous emission in the frequency interval 


v to v+dv is 
167 23 h c\2 dv 
-= —— 43Z?2 
de~ 57507 (ze) (5) F ) 


where « is the fine-structure constant; Ze is the nuclear charge; v is the initial 
velocity of the electron undergoing the transition; and h, m, and c have their 
usual meanings. It has become customary to write the quantum-mechanical 
cross-sections with which we shall be concerned in terms of Kramers’ formula 

dp =g dbx (2) 
where g is the so-called Gaunt factor. In cases of astrophysical interest the 
Gaunt factor is often close to unity, and it is frequently taken as such for all 
wave-lengths in approximate work. However, for long wave-lengths this is not 
a good approximation, and a more careful evaluation of g is required. 

We shall always work in terms of the spontaneous emission rather than the 
absorption cross-section. It is not difficult to obtain the latter from the former 
by applying detailed balancing in thermodynamic equilibrium, and it is in fact 
immaterial which we use if it is the Gaunt factor to which we always refer. 

Derivations of formulae for the emission cross-section for free-free transitions 
can be made using the usual methods of the quantum theory of radiation (3). 
The cross-section in the electric dipole approximation for a particle of charge Z,e 
to make a transition in which its velocity is changed from v;, to v; in the presence 
of an infinitely heavy particle of charge Z,¢ is 


ap= (5) azeze (72) | Kedar aiky|-da] S 


The direction of v, lies within the solid angle dQ, and the propagation vectors k; 


and k, are defined by the relation 
k=mv/h. 

From the classical equation of motion 

md?r/dt? = (Z,Z,e?/r*)r (4) 
one may obtain the relation 

Cllr Ya4(8, 8) = — TEES er *¥4(8 8) 

The matrix element on the right-hand side of this equation has been extensively 
discussed recently in connection with,the excitation of atomic nuclei by electro- 


magnetic interaction with charged particles (4), and, following these accounts, 
we may introduce the dimensionless function 





Fea 8)= EELS | Keylr*¥i4(0, 8k) dO 


The parameters n; and € are defined by the equations 
n= %Z,Z,(c/v;), f=17— 
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It is readily shown from the equations given above that 


= 9V3 Fn, £) (5) 


320% 





where, for the case of electrons (whose charge is to be taken as —e) in the field 
of a nuclear charge Ze, 
The remainder of this paper is concerned with the evaluation of the function 
on the right-hand side of equation (5). 


1. Formulae for the Gaunt factor 

1.1. The Born and Born-Elwert approximations.—The most elementary 
calculation of the Gaunt factor is that in which Born’s approximation is used. 
This is valid if the incident velocity is very large (|n;|<1) or if the energy loss 
is very small (i.e. pX%1, where p is the velocity ratio v,/v,). We set 


|k) =exp (ik.r)) 
and it is then easy to show that 
(kyr? ¥,"(8, 6)[k;) = a m7 Y,“(0,®) (6) 


where @, ® are the polar angles of (k;—k;). When this expression is inserted 
in equation (5) an elementary integration yields the formula 


gpa XV3 Reths\ _ V3), (e+! 
5 In n (py %) 7 in(2=*). (7) 


The conditions for the validity of Born’s approximation are not always 
realized under astrophysical conditions. Nevertheless the expression has the 
correct form in the infra-red limit p>1, which is an important region of application 
for many purposes. Elwert’s modification has a somewhat wider range of 
application since it takes some account of the modification of the wave-functions 
in a Coulomb field (9). The formula 

g® =". exp (27;)—1 (8) 
exp (277) —1 
is valid when either |é| <1 or |y,|<1, the restriction on 7; being much less stringent 
than in the Born approximation. 

1.2. Sommerfeld’s formula——The Born and Born-Elwert formulae are only 
valid over a restricted range of electron velocities, and outside that range it is 
necessary to use the correct Coulomb wave-functions in place of the plane waves 
of Born’s approximation. These wave-functions have the form 


|k;) =e722"%T'(1 +in;) exp (k;.r)®[ —in;, 1; (kyr —k;-7)]) 
Ie, = e224 (1 — in) exp (iky-F)Ling, 1; —i(hyr +ky-r)]) 
where the former contains only diverging and the latter only converging spherical 


waves superimposed on the plane wave. The functions ® are confluent 
hypergeometric functions. 

The integrals involving these functions have been calculated by Sommerfeld 
(6). He wrote each confluent hypergeometric function as a contour integral 
and then inverted the order of spatial and contour integration. In this way he 
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obtained the remarkably compact result 


- NAS p2nnj(e2xng__ 1)-1(e2%j — uiaf@f_ id Se eo 
g V3-"@e (e747 — 1)—1(e?"i — 1) lat x |F¢ inj, ings ts) } dO 


= 1/3. e? (e271 — 1) “1 (e247 — 1) ey |F(—in;, —ing; 1; Xo) |? (9) 
0 
where we define 
x=xysin?}0,  x%)= —4nin,/&, 


© being the angle of deflection of the electron. The function F is a hypergeometric 
function. 

Formulae equivalent to that of equation (g) have also been derived by 
expanding the wave-functions in spherical waves and summing over the partial 
cross-sections (7, 8). 

1.3. The formula of Menzel and Pekeris.—A completely different approach 
is that of Menzel and Pekeris (5). ‘Their method was to generalize by analytic 
continuation a formula given originally by McLean (10) for the Gaunt factor 
in line transitions averaged over a complete shell of principal quantum number n 
in a hydrogen-like atom. Menzel and Pekeris noted that, apart from certain 
factors in the normalization, the solutions of the Schrédinger equation for the 
discrete bound states and for states in the positive energy continuum can be 
written in the same functional form as a function of rand. The energy relation 
for a discrete level 





W = — «?Z?mc?/2n? 


can be inverted to give m as a function of W. If Wis positive, we have, on replacing 
it by 4mv?, where v is the velocity at a large distance from the atom, that 


n=1aZ(c/v) 


which is clearly —i times the parameter 7 used inthis paper. Thus, provided due 
allowance is made in the normalization for the fact that one can talk only about 
states within a given energy range in the continuum, rather than about a single 
state of specified energy, it is possible to derive formulae relating to bound-free 
and free-free transitions by expressing all energy dependence in terms of m in 
the formulae for bound-bound transitions, and substituting —7 for n wherever 
this is appropriate. 

In this way Menzel and Pekeris generalized McLean’s result for line transitions, 
and gave, in particulat, the expression 


gh? = my/g.ebmi(ePne— 1) He — 1 )-Ininle| 
* | [F(— in: + 1, ings 15 Xo)? — [F(t — ing +15 15 Mo) FP] (10) 
for free-free transitions. 
It appears that the direct equivalence of the Menzel and Pekeris formula 
with that of Sommerfeld has never been demonstrated. That the two expressions 


are in fact equivalent may be seen from the following considerations. We start 
from Sommerfeld’s formula (9). Consider the two relations (11) 


x4 F(ab; c; x)=a[F(a+1,b;c; x)—F(a,b; c; x)], 


(b—a)F(a,b; c; x)+aF(a+1,b; c; x)—bF(a,b+1; c¢; x)=0. 
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These may be used to transform the derivative occurring in equation (9) as 
follows 

d 


: : 3 d 5 : ae 
X95 |F(—im, — ings 1; Xo) [? =X» =— [F(—in;, —ings 13 Xo) F (ins nz; 15 %)] 
Xo Xo 
= inl F(—inj, — ing; 13 Xo) F (in + 1, ing; 13 Xo) 


— F(—ing+ 1, — ings 15 Xo) F (ini, ing; 15 X)] 
-i- [Fini + 1,275 15 Xo) F(—iny — Mme +15 15 Xo) 


— F(—in; +1, —ing; 15 %)F (ini, tnet+ 15 15 X)]. 
By means of the Kummer relation 
Y= F(— in; +1, — ing; 15 Xp) = (1 — Xp) PF (ini ing + 15 15 Xo) 
and its complex conjugate, the expression may be put in the form 


BOL {(x — Hq) H+ ED — [1 — mg) Hut EAP) 


where the asterisk denotes complex conjugation. Since the Gaunt factor is 
essentially positive, it is clear from equation (g) that the transformed expression 
that we have just derived is also real and positive. We may therefore put it equal 
to its modulus and write the Gaunt factor in the form 
g=my/3.emi(e?i— 1) (24 — 1) giny 
x lé|1( aa %y) Kut ng Xp2 [(1 — Xq) ut PP?) * |. 
Since (1 — 9) is real and positive, the factor 
(1 — x9 )euty) 
has unit modulus. Thus 
|(1 — x9) Hi tH? — [(1 — x9) Hit PP? ]* | = [P2 — [(1 — x9) 2 PP] I. 
Since the Kummer relation gives 
(1 — x9) C+ MY = Flin; ng +15; 15 Xo) 
we see that the Gaunt factor is given by 
ga my/3 .ermi(ermri— 1)-1(e"Y— 1) inglE| 
x | [F(— init 1, — ings 15 Xo) — [F(t — ing +15 15 Xo) Pb 
which is the form given by Menzel and Pekeris. 

1.4. The classical approximation.—A useful, and less complicated, formula 
for the Gaunt factor using Coulomb field wave-functions can be derived in the 
so-called ‘‘classical’’ approximation (4). ‘The W.K.B. approximation is used 
to solve the Schrédinger equation for the wave-functions, which can then be 
employed to evaluate the necessary matrix elements. In the classical limit of 
slow electrons the Gaunt factor becomes 


go=- os — "a5 Kya) } ™ (11) 


27 z=|é 





where K denotes a modified Bessel function. It is clear that this relation must 

be a confluent form of Sommerfeld’s formula. 
The same expression for the Gaunt factor can be obtained by studying the 
radiation emitted by a particle moving in a classical hyperbolic orbit in a central 








XUM 








246 I. P. Grant, Calculation of Gaunt factors Vol. 118 


Coulomb field (4,12). The agreement between the two approaches is in 
accordance with the Correspondence Principle, as we might expect in view of 
the fact that the W.K.B. approximation is essentially an expansion in powers of 
Planck’s constant. 

Finally, we shall take note of an important symmetry relation between the 
Gaunt factors in the two cases in which the particles involved have charges whose 
signs are the same or opposite. When the two charges have the same sign, the 
parameters 7; and 7; are both positive, and we may write 


(nine) =8(Inib Iny))- (12a) 
When they have opposite signs, the parameters ; and 7, are both negative, and 
8(ni 74) =8(Inil nel) exp (27 6)). (125) 


It is easy to verify these relations for all three formulae (9), (10) and (11). They 
have been used in the following section. 


2. Numerical values and the use of simple approximations 

The most extensive table that has been published so far is that of Thaler, 
Goldstein, McHale and Biedenharn (7) who have tabulated for like charges the 
function a, defined by the equation 


_2V3 4 

ae “— ry . (13) 
The greater part of Table I, which refers to the case of unlike charges, has been 
derived from the table of reference (7) by making use of the relations (12a, 5). 
Table I is divided for convenience into alternate regions printed in italic and in 
ordinary type, of which the three upper regions have been computed by this 
method. The accuracy quoted in (7) is to within a few units in the last place 
of decimals, and the same will be true of our table. 


TABLE I 


Electric Dipole Gaunt factors 


aa I‘o!I 1°025 1°05 IIo EES 1‘20 
0°OOI 2°885 2°419 2:047 1‘679 1-468 I +323 
oO'ol 2°906 2°420 2°049 1684 I°474 +330 
ov! 2°912, 2°430 2°072 ag 1°527 1°395 
0°25 2892 2°423 2082 1°760 1°583 1°469 
O's 2°815 2°365 2°046 1°755 1°604 1°507 
0°75 2°715 2°280 1°978 1°709 1°573 1°488 
1‘ 2°615 2°194 1906 1°654 1°530 1°452 
2°5 2°231 1868 1°629 1°440 1°352 1°298 
5°0 1°955 1°635 1°445 I-303 1+238 I-201 
7°5 1°807 1°518 I°357 1+239 1+187 I‘I57 
10°O t°707 1°446 I+303 I-201 I-I57 I°I3r 
12°5 1°634 1°394 1-266 I°I75 I+ 36 I‘It4 
15‘0 1°579 ¥°367 1-239 I°I57 I-I2r I‘r02 
20°0 1°497 1*302 I-201 I°I3I I-ror I:085 
25°0 1°442 1°265 I°I7§ I°-It3 T:088 I:073 
30°0 1-400 I+237 I-156 I‘ror r:078 I:065 
35°0 1-367 I°3I7 I-142 I-0g1r I:o71 1-059 
40°0 I°342 I-200 I-130 1-084 1-065 I:053 
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The fourth region, in the bottom right of the table, has been computed from 
2 a table constructed for like charges by Alder and Winther (13) using the 

formula (11). ‘Table II displays corresponding values of g° for unlike charges, 
h and it is readily ascertained that the classical approximation agrees with 
wh Sommerfeld’s formula to better than 0-4 per cent in the two lower regions of 
he Table I. Table II contains a further column in which g© has been calculated 


from an asymptotic expansion of equation (11) derived by application of the 
method of steepest descents. ‘This expansion is 
Cx oe | ae —4/3 
go~1+0°21775|€| 0:01312|é|-43 +... (14) 
valid when |€|>1. This formula gives remarkably good agreement with the 
b) exact value of g© even for |€| as low as 1. The first two terms of this expansion 
are given in reference (4) with somewhat fewer decimal places. 


"y ; : 
An asymptotic formula of the same character has been derived from 
equation (10) by Menzel and Pekeris. It is 
gMP a t +O°172¢ (1 +p-2)(p-} +p~?) Shs 5 i 
I, 
. 8 “— 
. —0-0496(1+ Fp +p-y(pt+p=yefae+... (15) 
15 


valid when |y,|>1 and p>1. This formula has not the same accuracy as that 
3) of equation (14), which is usable over a considerably wider range of photon 
frequency and electron velocity. 

The Born—Elwert approximation of equation (8) has been used to compute 
). Table III. It agrees with the exact numerical values to better than 0-2 per cent 


“i in the upper italic region of Table I. The Born approximation, if unmodified, 
- is satisfactory to this order of accuracy only for the lowest value of |nj, 
e 


\n;|=0-001, given in the table. 
In spite of this restriction, the Born approximation, modified or unmodified, 
has the same behaviour as all the other formulae in the infra-red limit p—>1, 


TABLE I (continued) 





1°25 1°30 1°35 1°4 1°6 1°8 

mt S. 

o'oo! I°212 I-I24 I‘O51 0-989 0-810 0:693 
o'ol I*220 I+133 I:06r I-000 0-824 0:708 
ov! I+296 1-218 I°I55 I‘I02 0951 0-856 
0°25 1°385 1°320 1‘268 1°226 i oe I*045 
O'5 1°438 1°387 1°346 1°314 1°230 1185 
0°75 1°428 1°383 1°348 1°320 1°249 1°208 
are) 1°398 E°357 1°326 1°303 5237 1°202 
2°5 1°261 1°235 I°2I5 I-198 1158 1-136 
5'0 I°I75 I‘157 I-I44 I+133 I:106 I-0gI 
7:5 I+137 I°T23 r-ir2 I*I04 1-082 1:065 
10'O T*Il5 I‘I03 IT-0904 1-087 Ir-069 1°054 
12°5 1-100 1-089 1-081 I:075 1°056 1°046 
15'0 r-088 I:079 I:072 1-067 1°049 T'042 
20°0 I:073 1-066 1-060 T+055 1042 1°034 
25°0 1-064 1-057 1-052 1°046 1°035 1°030 
30°0 I‘057 I-051 1°045 1°042 1032 1°'026 
35°0 I-O51 1°045 1°O41 1'037 1°029 1°024 
40°0 1-047 1°042 1°037 1°034 1°026 1°022 
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TABLE II 


Electric Dipole Gaunt factors. Classical approximation 


-&€ g°(Eq.11) = g© (Eq. 14) 
o'o oe) 

ov! 1°708 

o'2 1°502 

o°3 1°396 

O*4 1°344 

O's 1°306 

0'6 1°276 

0°7 1°252 

o's 1°232 

o'9 E255 

1'O 1‘201 1°205 
1°2 1°180 1183 
1°4 1°164 1°166 
1°6 I‘151 I°152 
1°8 1*140 I°I4I 
2'0 ¥*¥32 272 
4°0 1'084 1°0843 
6-0 1°0647 
8-0 10536 
10°O 1°0463 
15°0 1°0354 
20°0 1°0294 


Tas_e III 


Electric Dipole Gaunt factors. Born—Elwert approximation 





p "=e ge Ores oo VE O45 Ue Hors 

1‘O! 2924 2°925 2°932 

1°025 2°423 2°425 2°439 

1'OS 2°047 2°O51 2°076 

ry 1°679 1°684 1°726 

I'l5 1°468 1°475 1°S32 

52 1°322 1°330 1°398 

1°25 1231 1°222 1°298 

1°3 1°23 1°134 1°220 1°329 

£35 1°050 1°061 1°156 1°276 

I°4 0'9878 1000 1°103 1°232 

1°6 0°8084 0°824 0°952 1r°rts 1°246 
1°8 0°6907 0°708 0°856 1'047 1°194 
2°0 0°6057 0°625 0°790 1°003 1°161 
3°0 0°3822 0°*407 0°631 O°917 1°097 
4°0 0°2816 0°309 0°572 0°894 1°078 
5°0 0°2235 0°252 0°545 0886 1:069 
6°0 o'1855 0216 0°532 0882 1°065 
8-0 0°1385 o'171 0°520 0878 1°060 
10°0 o'1106 0°144 0°517 0°876 1°058 
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that is to say 


; > 
A (16) 
7 a= Tt 


on 
» 


The singularity of the Gaunt factor in the limit p->1 is connected with the infinite 
scattering cross-section of the Coulomb field. It is to be expected that this 
singularity will be strongly modified in atoms other than hydrogen by screening 
of the atomic electrons. We shall now consider this modification. 


3. Screening in the infra-red limit 

The fact that the Born approximation has the correct behaviour in the infra-red 
limit suggests that the method be applied to calculate the Gaunt factor in this case. 
The cross-section for free—free emission in the non-relativistic dipole approxima- 
tion is given by equation (3). Suppose that the screened potential function is 
written in the form eZ(r)/r, so that Z(r) is the effective charge for the potential 
at radius r. Then from the equation of motion one obtains 











e Z(r 
(kr ¥,"(4, $)|k;) = — atm (Ks 5 ( ) ¥,"(0,$)k) 
which leads to the formula 
kik; i d {Z(r) > alo 2 
ee. oe mie E u AY Ik. Q2. 7 
5 G@nyry3 >) (ky = ( r ) Y, ( rp) ) d (17) 








In the screened field the analogue of equation (6) is 
d (Z(r) : _ me eae Z(r) 
—{—}. Y,“(0,¢)] k;) =4m1 Y,"(Q Kr 
5( r "O9) Y site ), (A 5(4 )) ad 


which, after integration by parts, equals 


f 








~ 4miY,*(Q’) [ Z(r) sin Kr. dr 
J0 


where j,(AKr) is a spherical Bessel function of order unity, K=k;—k,, and Q’ 
denotes the polar angles of K. ‘Thus 


” Z(r) 
|. 7 sin Kr. dr 


4 





g= Shi ik | “dQ. 








Since K2=k?+k,?- site and dQ=2z7dy, this result can be written 
> V3p ty | I. | - ") 
ve 0 


This expression is equivalent to, al somewhat different from, the expressions 
discussed by Bethe and Heitler (19) and Kirkpatrick and Wiedmann (20). 
These involve the substitution of Z—F(q) for Z in the formula for the cross- 
section before integrating over the angles, where 


sin Kr. dr| i (18) 








F(q)= | plrye e'4-* dr 


is the atomic form-factor, and fq is the recoil momentum of the atom defined 
by the relation 


q=k,—k,;-k=K-k, 


in which /K is the photon momentum, so that k= 2z/c. 
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The expression (18) may be evaluated either numerically as it stands, or else 
by means of an analytical approximation. Such an approximation is 
Z(r)= E a,e-F0" 
where Z(0)=Z, and 
lim Z(r)=Z* 
where Z* is the degree of ionization of the atom. The Gaunt factor has the 
general form 


= V3 a2 oe + (k; +k)? _ 4B,,*kiky 
ant 2 [l* B+ (Ry B+ RFU + ek) 


an Batt (hth)? _p ay Butt (hi+hy) 
+2 5 gttea| Btn ep By Pm nas male (9) 


In the case of an ion, the contribution of the unscreened part of the field may 
be written separately, and if we put a=Z*, By=0, 


e-(Z) eae (Z) Sara oe 


where g’ is of the form of equation (19). If all the coefficients in the formula 
for the potential are zero save that for m=o, the result is just the formula for 
Born’s approximation for free-free transitions in the field of a bare charge Z*. 

Although this paper is concerned mainly with positive ions, the theory of 
screening effects in this section can be applied to neutral atoms also, provided 
that an effective potential can be given in which the electron is considered to 
move. ‘The potential function is obtainable from, say, the Thomas—Fermi 
theory of the atom or else from a Hartree self-consistent field calculation. In the 
latter case, the required function Z(r) is one-half of that usually denoted as 
“Total 2Z,,”’ in the literature (18). 











4. Average Gaunt factors 

It is often necessary in astrophysical problems to compute a cross-section 
which is averaged over a statistical distribution of electron velocities. This has 
been done, for example, by Berger (14) for a dilute hydrogen plasma, and also 
by Menzel and Pekeris (5). In this section we shall discuss a simple method of 
making such calculations for the Maxwell distribution. 

To begin with, let f(v)dv be the probability of finding an electron in the 
velocity range v tov+dv. Then from the definitions (1) and (2), the cross-section 
d®(v) for emission of a photon in the frequency range v to v+dy, is given by 


167 350/ h \* dv f{ (* awhe¥ roo ae 
ao(v)= 25 002(£) 24 |" een (E) aseddn} = {fafean} (ax) 
where vp is the minimum velocity for which the frequency v may be emitted. 


The average Gaunt factor G(v) is defined as the ratio of formula (21) to the 
corresponding expression obtained from Kramers formula, so that 


Gr)= |" eoonnenS |" fay (22) 
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If we now suppose that f(v) is a Maxwellian distribution of velocities we may 
write this in the form 
G(v)= [- g(r, eierme™” odo, | | emmeat o,do;= |" guy x)e--™ de (23) 
Ve vr) u 
where 
u=hv/kT, x=mv?/2kT. 
The absolute temperature is denoted by T and k is Boltzmann’s constant. An 
alternative expression to (23) is obtained by writing 
y=x—-u 
so that y is just the final energy of the electron divided by kT. In this case 


Glw)= |“ soyeray (24) 


where g(y) is to be evaluated at the values of p and 7; given by 


my Per u\ "2 a pet 12 
“in ind yru 


The integral (24) can be evaluated numerically by any suitable method of 
quadrature. One such method is given by Chandrasekhar (15); it is based on 
the evaluation of the integrand at the zeros of the Laguerre polynomials, the 
values thus obtained being added with appropriate weights to obtain the value 
of the integral (see Table IV). This method is, to anticipate the examples of 
the next section, extremely rapid in application. Only three points were needed 
to obtain three significant figures in most cases. 


TABLE IV 


Quadrature for the mean Gaunt factor 


«oo ™m 
Ge)= | e-veody= Yanelon) 
. I 
m=I y=1 a,=1 
m=2 ¥1=0°5858 a, =0°8536 
Y2=3°4142 a,=0°1464 
m=3 ¥1=0°4158 a,=oO°71II 
Y2=2'2943 A, =0°2785 
¥3 = 6-2899 a3 =0°O104 
m=4 V1 =0°3225 a, =0°6032 
Yo=1°7458 dy =0°3574 
¥3=4°5366 a3 =0°0389 
¥a=9°3951 @4=0°0005 


It is, of course, possible to obtain an analytical formula for G(v) from the 
asymptotic formula (15), as Menzel and Pekeris have done, by expanding it 


in the form 
hy 1/3 2y 
gxi+or728( 5) (1424...) +... 


2m*mZ"e* 
for y <u, a condition equivalent to taking p>1. The integration over y may then 
be performed explicitly to give 


34, 15 \ 
Gv) 1401728 (a) (142+...) +.... 


2n?*mZe! 
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It is clear from the method of derivation that this is valid only if u is so large that 
the second term in the bracket is very small in comparison with the first. This 
imposes a severe restriction on the use of the result. 

A similar but more complex expression may be derived from the asymptotic 
formula (14). This suffers from the same disadvantages. The use of both results 
is therefore somewhat restricted, and the advantage of having an explicit formula 
is not sufficient to recommend them for practical calculations as a general rule. 
The numerical method described above is usually much to be preferred. 


5. Some numerical examples 


The object of this section is to illustrate by means of a few examples the way 
in which the methods suggested in this paper may be used to obtain Gaunt 
factors both for particular values of the photon frequency and electron velocity, 
and also for average values taken over a Maxwellian distribution. They have 
been chosen more with the idea of demonstrating the ease with which such 
calculations can be performed rather than with a particular physical situation 
in mind, and in the case of the examples involving screening, the data immediately 
at hand determined the choice that was made. 

Example 1. Hydrogen Z=1. 

A=c/v=500A, T=o'5eV, u=49°6. 
The mean Gaunt factor in this example is given in the table of the paper by 
Berger (14) already cited. His value is G(v)=1-10. 


We shall use the numerical method of the last section, the weights being 
obtained from Table IV. 


One point formula G(v)=g(1-0) = 1-079. 
Two point formula G(v)=0°854 g(0°586) + 0-146 2(3°414), 
8(0°586)=1-065, = -8(3°414) = 1-131, 
G(v)= 1-074. 


Three point formula G(v)=0°711 g(0°416) + 0°279 g(2°294) + 0-010 g(6-290), 
&(0°416)=1-057, —_g(2294)=1-110, (6-290) = 1-173, 
G(v)=1°073. 
In all the values of y occurring in this example |é| was larger than unity and this 
enabled the use of the asymptotic expansion of the classical approximation (14) 
in the calculation of g(y). 

The corresponding calculation by Berger involved the computation of g for 
two hundred incident electron energies. These values of g were integrated 
over a Maxwell distribution using Simpson’s rule. The comparative simplicity 
and the ease of our method are evident. 


Example 2. Hydrogen 7=1. 
A=10'A, T=100eV, u=0°0124. 
In this example, taken from the other end of the range covered by Berger, we 
have 








~ * 


0:0124\1/2 00062 
p={1I+ ; wit 


Clearly p always remains close to unity, corresponding to the infra-red limit in 
which the Born approximation is applicable, 
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One point formula G(v) =g(1-0) = 3°184. 

Two point formula __g(0-586) = 2-894, 2(3°414) = 3°860, 
G(v)=3°035. 

Three point formula g(0-416)=2-701, g(2°294)=3°643, (6-290) =4°198, 
G(v)=2-979. 


The rate of convergence with increasing number of points is not particularly 
rapid in this case because of the influence of the logarithmic singularity which 
makes g(y) much more strongly dependent on y than in Example 1. It is possible 
to obtain a better estimate of G(v) by the well-known method of exponential 
extrapolation (21), that is to say by taking 


G(v) = Ga(v) —{ Ga(v) — Ga(v)}?/{Ga(v) — 2G4(v) + G,(v)} 


where G,,(v) denotes the value of G(v) given by the m-point integration formula. 
With the values obtained above one finds 


G(v)=2°945 
for the present example. The difference between our value and the value of 
2°817 given by Berger is almost certainly due to our use of the Born approximation 
in a region in which he claims his numerical method is most accurate. 

These two examples illustrate how simple the calculation of mean Gaunt 
factors is in two extreme cases in which the classical approximation and the 
Born approximation can be used. Intermediate cases can be treated by 
interpolation in Table I, linear interpolation in p and logy, ( — n;) being satisfactory. 

We now turn to the estimation of the effect of screening on the Gaunt factor, 
following the method of Section 3. For this one may use any suitable potential 
function Z(r), and two such functions will be discussed. The temperature in 
both cases will be 1oo eV, and the wave-length 1o* A, as in Example 2. 

Example 3. Complete screening. 

We shall first use the universal Thomas—Fermi potential function for the 
field of an isolated atom as tabulated by Gombas (16), approximating to it by 
a series of exponential terms. The constants are determined by the method of 
Prony (17). 

Suppose we wish to represent Z(r) by, for example, a three-term series of 
exponentials 

Z(r) = ape—F" + x,e—Fi? + aye Far, (25) 


Let Z,,=Z(r)+nh), where n may take the values 0, 1,...5. Ifthe representation 
of Z(r) were exact, one could find constants A, B, C such that 


AZ,+ BZ,+CZ,+Z)=0 | 
AZ, + BZ,+CZ,+Z,=0 7 (26) 
AZ, + BZ,+CZ,+Z,=0 | 
Clearly if e~*", e-*", e-*" are roots of 
A#®+ B’e+Ct+1=0 (27) 
the equations (26) are satisfied identically whatever the value of 7». Thus we 
solve (26) for A, B, C with the known values of Z,, and obtain the 8; by solving (27). 


The «; may then be determined either by fitting (25) at three points with the 
values of 8; already obtained or else by the method of least squares. 
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Two exponentials were used to fit the Thomas—Fermi potential, the parameters 
obtained being 
A = 0°309Z, By = 0°313249/ap, 
% =0°691Z, B, = 1°316Z1/ap, 


where dy is the first Bohr radius in hydrogen. For large r the Thomas—Fermi 
potential function Z(r) behaves as r~*, so that the Prony series decays more 
rapidly. In this example, the deficiency, which increases with increasing r, 
has only reached to per cent at a point at which Z(r) has decreased to 3 per cent 
of its maximum value. The corresponding error in the radial integral is small 
since Z(r) is well represented in the important part of the range of integration. 

Using a three-point quadrature formula we obtain the results displayed in 
Table V. 


TABLE V 
Comparison between screened and unscreened Gaunt factors 
a(y) 
y Screened (Thomas—Fermi field) Unscreened 
Z=4 Z=82 Any Z 
0416 0°334 0063 2°701 
2°294 0°726 0°243 3°643 
6-290 0991 0°437 4°198 
G(v) 0°450 O°117 2°979 
Example 4. Neutral Be. Z=4. 
A=10'A, T= 100 eV, u=0°0124. 


It is of some interest to compare the Gaunt factors computed from the 
Thomas—Fermi potential with those computed using self-consistent fields. 
In the case of neutral beryllium the self-consistent field is that computed by 
Hartree and Hartree (18). The representation of their values used in this 
example makes use of the fact that for large r the screening due to the shell of 
largest radius has the form Z(r)~exp(—+/(2e).7), where « is the eigenvalue 
in atomic units. For the case of neutral beryllium in its ground state the 
configuration is (1s)?(2s)*, and the 2s eigenvalue has the value 0-5764. The 
asymptotic behaviour of Z(r) is therefore exp (— 1-07377r), and the approximation 
used in this example to represent the self-consistent field proceeds by subtracting 
this from the numerical values, and using Prony’s method with two exponentials 
for the remainder. ‘The two exponents yielded by this process are found to be 
complex, and the approximation obtained has the form 


Z(r)/Z =e-1 073" — 9-557 1e—2 352 sin (230427). 


The two conjugate complex terms have been combined in the last term of this 
result. 

It is best to work with real expressions for g, and for this reason the same 
argument as leads to equation (19) has been applied to the expression 


Z(r) =ae-®r — Ae—¥r sin Cr. 
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The resulting formula is 














g= ae | {2 ™ B? + (kh; + ky)? ms 4PB°kik, ] 
anZ? BP+(k;—Ry)? {BP + (ky + hy} B? + (ki —hy)?} 
a 2%82ABC In (k; + k,)* + 2(B*— C?)(k, + k,)? + (B2 + C2)? 
(6? — B? + C?)? + 4B?C? (k;— ky)* + 2(B? — C?)(k;—k,)? + (B2+C?) 
(Seema) ae tcc b+ cop 
BP + (Rk; + ky)? . (A; +k,)* + 2(B? — C?)(k; +k,)? + (B2+C?)P 





7 (B?— C2)(k,— ky)? + (B® + C2)? : Be ~C?) 
(k,—k,)' + 2(B*— C*)(k,—kyP+(BP+O| | 4BC 
f f 
AAPOR PEON) 5 ans eth EBC 
4B°C? + (B?— B?+ C2 2BC 
_, i—ky)P + BP a 
2BC ; 


The numerical values obtained from this formula for the neutral Be self-consistent 
field are given in Table VI. 


+20. 








—an 





TABLE VI 


Screened Gaunt factors for neutral beryllium using the 
self-consistent field of Hartree and Hartree (18) 


y a(y) 
0°416 0°261 
2°294 0°546 
6-290 o°751 
G(v) 0°345 


As is to be expected, the value of G(v) in Table VI computed from the 
self-consistent field is smaller than that obtained from the Thomas—Fermi field, 
G(v)=0-450. The profile of Z(r) in the latter case lies consistently above that 
for the self-consistent field as is well known. The difference between the two 
values of G(v) is in fact insignificant by comparison with the remarkably large 
reduction from the unscreened value of the Gaunt factor, G(v) = 2-979. 

Example 5. Bet+. Z=4. 

A=10'A, T = 100 eV, u=0°0124. 

It is instructive to make a further comparison of the last example with the 
case of doubly ionized beryllium, for which the self-consistent field is again 
given by Hartree and Hartree (18). The nuclear charge is screened in this 
case only by two 1s electrons, and a sufficiently good approximation to Z(r) 
appears to be 

Z(r)=2[1 + exp (— 4°7637)]. 
This approximation is correct at r=o and for r greater than about 1-0, and is 
nowhere in error by more than about 5 per cent. 

We write equation (20) in the form 


&(y)=82(y) + 8:() +8'(y) 
where g,(y) denotes the term containing a factor (Z*/Z)?, g,(y) the term containing 
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(Z*/Z), and g’(y), has the same meaning as was assigned to it in equation (20). 
With a three-point quadrature formula the results are as given in Table VII. 


TABLE VII 
Screened Gaunt factors for Be++ using the self-consistent field of Hartree and Hartree (18) 
y go(y) £9) g(y) o(y) 
0416 0676 0060 0006 0°742 
2°294 o'9gI1 0-190 0°044 1°145 
6-290 1°050 0°305 o°0gI 1°446 
G(v) 0°745 0°099 0°O17 0861 


The Gaunt factor in Table VII is over twice the value of that for the neutral 
atom, but is still a factor of nearly 4 down on the unscreened value, G(v) = 2-979. 
The value of the Gaunt factor in fact differs by only 10 per cent from the value 
required to make the cross-section equal to that for a hydrogen-like system of 
charge Z* =2. 


Conclusions 


The use of Tables I, II and III of this paper, combined with the available 
asymptotic formulae, enable one to make for hydrogen-like systems a rapid 
calculation of the monochromatic Gaunt factor for given electron velocities. 
The rapid methods of Section 4 enable one to verify the tables given by Berger (14) 
for the mean Gaunt factor in a dilute hydrogen atmosphere within tolerable 
limits over the whole of his range of temperature and wave-length. This range 
may be easily extended using the same technique, a useful feature in view of the 
fact that Berger’s tabular values of temperature and frequency have to be scaled 
by a factor Z? to give Gaunt factors for a hydrogen-like atom of charge Z. 

The estimates of the effect of screening that have been given in this paper 
demonstrate how important it is to take it into account when dealing with 
free-free transitions in real atoms and ions. The Thomas—Fermi field gives 
a useful means of estimating the cross-section for neutral atoms when no 
self-consistent field is available, and a first approximation to the cross-section 
for a positive ion can be obtained by considering it as if it were a hydrogen-like 
system with a charge equal to the degree of ionization. This last approximation 
is, however, likely to be somewhat inaccurate if the system is only slightly ionized. 

Acknowledgments.—The writer would like to acknowledge the benefit of 
consultation with Dr D.W.N.Stibbs to whom he is indebted, in particular, 
for suggesting the use of Prony’s method in the calculation of the effect of 
screening on the Gaunt factors. 

Theoretical Physics Division, 

A.W.R.E., Aldermaston, 
Berkshire: 
1958 March 12. 








) 





7) 


No. 3, 1958 for free—free transitions near positive ions 257 


References 
(1) H. A. Kramers, Phil. Mag., 46, 836, 1923. 
(2) G. Wentzel, Z. Phys., 27, 257, 1924. 
(3) W. Heitler, The Quantum Theory of Radiation, p. 249, 1954 (Oxford). 
(4) K. Alder, A. Bohr, T. Huus, B. Mottelson, and A. Winther, Rev. Mod. Phys., 28, 
432, 1956. 
(5) D. H. Menzel and C. L. Pekeris, M.N., 96, 77, 1935. 
(6) A. Sommerfeld, Atombau und Spektrallinien II, 1951 (Vieweg, Braunschweig). 
(7) R. M. Thaler, M. Goldstein, J. L. McHale and L. C. Biedenharn, Phys. Rev., 102, 
1567, 1956. 
(8) L. C. Biedenharn, Phys. Rev., 102, 262, 1956. 
(9) G. Elwert, Ann. Phys., 34, 178, 1939. 
(10) L. McLean, Phil. Mag., 18, 845, 1934. 
(rx) Bateman Manuscript Project, Higher Transcendental Functions I, 103, 1953 (New 
York, McGraw Hill). 
(12) L. Landau and E. Lifshitz, Classical Theory of Fields, p. 200, 1951 (Cambridge, 
Addison-Wesley). 
(13) K. Alder and A. Winther, CERN report T/KA-AW-1, 1954. 
(14) J. M. Berger, Ap. }., 124, 550, 1956 ; Phys. Rev., 105, 35, 1957. 
(15) S. Chandrasekhar, Radiative Transfer, p.63, 1950 (Oxford). 
(16) P. Gombas, Die Statistische Theorie des Atoms, 1949 (Springer-Verlag). 
(17) E. T. Whittaker and G. Robinson, The Calculus of Observations, p. 369, 1948 (Blackie). 
(18) D. R. Hartree and W. Hartree, Proc. Roy. Soc., 149, 210, 1935. 
(19) H. A. Bethe and W. Heitler, Proc. Roy. Soc., 146, 83, 1934. 
(20) P. Kirkpatrick and L. Wiedmann, Phys. Rev., 67, 321, 1945. 
(21) D. R. Hartree, Numerical Analysis, p. 32, 1952 (Oxford). 








ATOMIC AND MOLECULAR NEGATIVE IONS IN STELLAR 
ATMOSPHERES 


L. M. Branscomb* and B. E. }. Pagel 


(Communicated by the Astronomer Royal) 
(Received 1958 March 28) 


Summary 


Late-type stars show continuous absorption from other sources besides 
the H~- ion. These sources begin to appear in the Sun, and increase in 
strength towards later spectral types, especially in carbon stars. On the basis 
of recent theoretical and laboratory work on atomic and molecular negative 
ions, the concentrations of C~-, O-, OH- and CN°~ are roughly calculated in 
giants and dwarfs of various compositions and their contributions to the 
continuous absorption coefficient assessed. O- is generally unimportant; C- 
is likely to be important in stars with very high C abundance such as R Cr B; 
OH~- may be important in the late M and S stars; and CN- perhaps 
contributes more than C; to the violet absorption in late N stars. Observa- 
tional verification of these conclusions will require accurate spectrophoto- 
metry over a wide range of wave-length, in addition to further laboratory 
work. 





1. Introduction.—Continuous absorption by the H~ ion, whose astronomical 
importance was first suggested by Wildt (1939), accounts for most of the 
observed features of the continuous spectra of the Sun and other stars of 
intermediate temperature. However, additional sources of continuous absorption 
appear to be present in the cooler layers (below 7000 °K) of the Sun (Michard 
1953), and in late-type stars (Keenan 1954) shortward of A4500A. A very 
heavy depression of the violet in late-type carbon stars shortward of A 4300 
has been attributed to continuous absorption by the C, molecule (Phillips and 
Brewer 1955; McKellar and Richardson 1955). 


Early attempts to explain these additional absorption continua by means of 
negative ions (Wildt and Chandrasekhar 1944) were handicapped by lack of the 
necessary laboratory data, and their failure has led to the suggestion that they 
should be interpreted primarily in terms of absorption by neutral molecules, 
either in continuous absorption to repulsive states or in the blending of the 
lines of many band systems (Wildt 1957). ‘The purpose of the present note is 
to evaluate the contribution of atomic ions other than H~ and to recall attention 
to certain further possible contributors, namely molecular negative ions, which 
were briefly mentioned as possibilities by Swings (1951). Recent laboratory 
studies of photodetachment spectra of both atomic and molecular negative ions 
make it possible to calculate, in at least a semi-quantitative manner, the 
circumstances in which the ions may make significant (if not dominant) 
contributions. 

The general features of photodetachment spectra are described by Massey 
(1950) and Branscomb (1957). The spectra of atomic and molecular ions are 
similar to the ionization continua of corresponding atoms or molecules, except 


* On leave from National Bureau of Standards, Washington 25, D.C. 
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that (a) the ‘‘ ionization potential ’’ of the negative ion (the vertical detachment 
energy in the case of molecular ions) is of the order of only a few electron volts ; 
(5) the continuum has no adjoining discrete line or band spectrum; and (c) the 
shape of the photodetachment spectrum near the threshold is markedly different 
from that of the ionization continua of neutral atoms or molecules. The most 
astrophysically interesting negative molecular ions are OH~-, in M and S stars, 
and CN~- and C,~, in carbon stars. The other molecules of relatively high 
stellar abundance either do not form negative ions at all (such as N,, CO, and 
NO) or are weakly bound and of relatively low abundance (such as O,~, CH- 
and possibly NH7-). 

The character of the wave-length dependence of the continuous absorption 
coefficient is determined by the electron configuration of the specific atomic 
or molecular ion. While photoionization cross-sections in general decrease 
initially from a finite value at the threshold, the photodetachment cross- 
section vanishes at the threshold. In the case of atomic negative ions, the 
cross-section rises approximately as vE*? for /=o and as vE@!-)? for />1, where 
E is the ejected electron kinetic energy (E=Av—E,), and / is the angular- 
momentum quantum number of the valence electron. Thus the H~ cross-section 
rises gradually from threshold (as £**), reaching a maximum at a photon energy 
roughly twice the electron affinity. The absorption of O- and C-, however, 
would rise as E1?, with initially infinite slope, giving an abrupt absorption edge. 
In both cases one expects a reasonably slowly varying cross-section, without 
pronounced structure. 

Somewhat analogous rules have been shown by Geltman (1958) to govern 
the behaviour of molecular negative ions near threshold. Here one expects the 
abrupt absorption varying as E? (except for consideration of the vibrational 
energy distribution) for heteronuclear molecules with the valence electron in a 
o or 7 orbital (such as CN~ and OH-), and for electrons of homonuclear 
molecules in those states with ungerade symmetry. Homonuclear molecules with 
valence electron orbitals of gerade symmetry (such as O,~ and C,~) will have 
photodetachment cross-sections rising as E**, ‘These general rules have been 
found to be in accord with the limited laboratory observations of OH~ and O,~. 

Laboratory spectra have been measured only for atomic ions of hydrogen 
(Branscomb and Smith 1955), oxygen (Smith and Branscomb 1955) and 
sulphur (Branscomb and Smith 1956); and for the molecular ions OH- 
(Branscomb 1957) and O,~ (Burch et al. 1958). However, with recent 
determinations of the electron affinities of the other ions of interest, and with 
the above rules governing the shapes of the cross-sections, one can make 
qualitative predictions about the other ions. ‘Table I gives the values of the 
photon energy and the optical wave-length which we shall take to be the threshold 
for the absorption coefficient of the several ions of interest. Several qualifying 
comments are required here. The binding energies given for H~ and O~ are 
accurately known, as are the photodetachment cross-sections. The value given for 
the C- threshold is also experimental, although not based solely on a photo- 
detachment experiment. Therefore, although the threshold energy and 
approximate shape of the cross-section are known, the absolute magnitude of the 
cross-section is not known. This. statement also applies to CN~ and C,~ in the 
table. The absolute cross-sections have been measured for OH~ and also for 
O,-, which does not appear in the table. 
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TABLE | 
Neutral species Negative ion 
Ground Ground Detachment Threshold Energy Ref. 
state Ou state QO. energy,eV wave-length dep.* 
H *Sij2 2 |; oe I 0754 16450A £E3/2 a,b 
O "Po 1,2 9 O- *Pij2,32 6  1°465 40-005 8460 Ele c,d 
Cc *Po ie 9 C- 4Sgi 4 112 +0°06 11070 Eve e 
OH “Ili 32 4 On” *2* ¥ £76 6970 Ell e 
CN *#Zt 2 CN- 15+ I 31+0°1 4000 Ev f 
Cc “I, or 3 C.- °%,+ - E3/2 
2 ‘E+ 1 2 g a 3900 g 


* E is the kinetic energy of the ejected electron: 

hc/A=EA (XY) +E. 

References: 

Branscomb and Smith 1955 

Chandrasekhar 1945 

Smith and Branscomb 1955 

Branscomb, Burch, Smith and Geltman 1958 

Branscomb 1957 

Bakulina and Ionov 1954 

g. Honig 1954 


mean oe 


There is one further complication. The electron affinity of a molecule is not 
equal to the vertical detachment energy of the molecular ion, if the formation 
of the ion leads to a shift in the equilibrium internuclear separation. We need 
the electron affinity in order to calculate the equilibrium concentration 
of the ions in stellar atmospheres. We need the vertical detachment energy 
in order to calculate the wave-length at which the absorption will set in. For 
the purposes of this paper we shall assume that the affinity and vertical detachment 
energy are identical to one another for OH and CN. In the case of OH 
there is experimental evidence suggesting that this is true (Branscomb 1957). 
However, even if it is not true, the affinity is larger than the value given in 
Table I, probably by one vibrational quantum of OH-, (i.e., EA (OH) ~2.2 eV), 
in which case the OH~ would be about 10 times more abundant than our 
calculations suggest. For CN- the value given is the affinity itself, and if 
this is not equal to the vertical detachment energy, only the wave-length at which 
the absorption occurs will be in error, and not the CN~ abundance. 

One must also keep in mind that laboratory photodetachment spectra 
probably are characteristic of absorption by ions without much vibrational 
excitation. In a stellar atmosphere, the higher vibrational states of the negative 
ion will be well populated, giving the absorption spectrum a rather different 
appearance. Thus for OH~-, for example, we expect the laboratory photo- 
detachment spectrum to be modified in two ways. First, the apparent threshold 
for absorption in the Av=o sequence may be shifted to longer wave-lengths. 
This will be true if w,” >w,’, i.e. if the vibrational energy levels of the negative 
ion are more widely spaced than those of the neutral molecule. Second, the 
higher temperatures in the star will give rise to a Av = — 1 sequence* which might 

* By a transition (v’, v’) in a photodetachment spectrum, we mean the continuous 


absorption of a negative ion XY~ (v”) in vibrational state v”, giving rise to a free electron and 
a neutral molecule XY (v’) in vibrational state v’. 
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at room temperature be absent, and which gives an additional continuous 
absorption threshold one vibrational quantum to the red of the main (Av=o) 
absorption edge. Because of such effects an absorption edge in a stellar spectrum 
might lie several hundred angstroms to the red of the laboratory absorption edge. 

2. Abundances of negative ions.—In the tables which follow we have 
attempted to estimate for each of these ions both the ratio of abundance of the 
negative ion to H-, and the ratio of its abundance to CN, taken as a representative 
neutral molecule of relatively high abundance and spectroscopic prominence. 
For carbon stars we have attempted to estimate the relative concentrations of 
CN~- and C, and C; molecules in the initial states appropriate for the absorption 
of the Swan bands and C, continuum, respectively. For simplicity in these 
necessarily qualitative calculations, we have considered several stellar models 
(Table II) from the calculations of de Jager and Neven (1957). Retaining their 
notation for ease of reference, we take compositions Ia, Ib, and Ic as approximately 
representative of normal and S stars, and take composition IIIb for carbon stars. 
For the relationship between spectral type, temperature in the reversing layer, 
electron pressure, and total gas pressure, the data given by Allen (1955) were used. 


TasBLe II 
Stellar composition models of de Jager and Neven 


Model log N (A) log N (C) log N (N) log N (O) Remarks 


Ia 00 —3°4 —3°4 —30 Fairly typical, but un- 
certain by factors of 
2 or more 

Ib oo —3°4 —3°4 —2°3 oxygen-rich 

Ic o"o —3°4 —3°4 —3°4 oxygen-poor 

IIIb oo —2°0 —2°'0 —3°0 carbon-rich 


Electronic statistical weights for the molecular negative ions (‘Table I) were 
found by assuming that the ground states of the negative ions have the same 
configuration and term as the ground states of iso-electronic molecules. (The 
energy dependence of the absorption cross-section near threshold rests on this 
same assumption.) The electronic statistical weight is 


Ou ints Gon . Oy . O,, ’ 


where ois is $ for homonuclear molecules, 1 for others; Q, is 1 for & states 
and 2 for others; and Q,, is the multiplicity for & states and other states with 
Hund’s case (6) coupling. Q,, is the number of values of Q in the state for 
Hund’s cases (a) and (c) (Gaydon 1947; Glasstone 1945). This Q,, is smaller 
by the factor (2A +1) than the value used by Pecker and Peuchot (1957), which 
was taken from the work of Hunaerts (1947) and Roach (1939). Although de 
Jager and Neven used the dissociation constants computed by Pecker and 
Peuchot, the magnitude of the error (a factor 3 for molecules with II ground 
states) is not sufficiently large to justify recomputing the equilibrium molecular 
compositions in view of the considerable uncertainties that remain in the 
abundances of the relevant elements. 

We make one further simplifying assumption. The rotational and 
vibrational constants of the molecular ions are assumed to be unchanged from 
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the parent neutral molecule, permitting us to use Saha’s equation in the form 
appropriate for atoms. This assumption is an approximate consequence of our 
assumption that the ground-state internuclear separations in the ion and neutral 
molecule are the same. In any case the actual molecular constants for the 
negative molecular ions are unknown. 

Since in stars the electron pressure is so low that the negative ion concentration 
will always be smaller than the concentration of the corresponding neutral 
compound, it is clear that the atomic ions will be important only in atmospheres 
of intermediate temperature if at all. In stars of types F to K the molecular 
concentration is sufficiently small that the continuous spectrum is certainly not 
determined by the blending of band spectra, so one need only calculate the 
ratio of the C~ and O~ absorptions to that of H~ to get a measure of their 
importance. ‘To the approximation required by the assumption of the preceding 
paragraph, this ratio is found from Saha’s equation: 


P(X)p(e) _ , Q(X) (27m\*? 02 50 BAO 
BY 2G (a) ATM rons, 

where Q,, is the appropriate electronic statistical weight; EA(X) is the electron 

affinity (in electron volts) of molecule X, and © is the decimal Boltzmann factor 

5040/T. ‘The ratio of the concentrations of two negative ions is then 


p(X-) _ ,, O(X-)O(Y) a | 
p(X) log O(X)O(Y>) + [EA(X)-EA(Y )]0 + logy) 

Since C~ and O~ are very unlikely to be of any importance unless the hydrogen 
content of the star is unusually low, the ratios N(C~)/N(H~) and N(O-)/N(H-) 
were computed for the most favourable cases in Table II. The ratio of the 
C~- to H~ concentration in dwarfs of composition Ic (low oxygen) and in giants 
of composition IIIb (carbon stars) is given in Table III, in which the ‘‘C ”’ 
classifications are taken from Keenan and Morgan (1941). The C~ is seen to 
be quite negligible in the main-sequence stars of normal composition. However, 
in the earlier giants the ratio of N(C~) to N(H~-) reaches approximately twice 
the fractional concentration of carbon in the atmosphere. Thus in G and K 
giants, which are too hot for molecular bands to dominate the spectrum, one 
may find cases of stars sufficiently deficient in hydrogen for C~ to contribute to 
the continuous opacity, provided that the C~ photodetachment cross-section is 
not abnormally small. R Coronae Borealis (Berman 1935) is an example of an 
object whose infra-red spectrum should be examined carefully for a break in 
the continuum near 1-1 microns. (See Section 4.2 below.) 











log 


TaB_e III 


N(C-)/N(4-) for composition Ic dwarfs and IIIb giants 


1) 2°00 «I'70—'—«=*dT'*4O”—é‘«‘~T'2O—~CT'1000—S_—s«O0'BKOD 
Camm. Se f spectral type Ms M3 Ks-Mo Kz2 Gs F4 
dwarfs ? log P 5°0 5°0 5"0 4°95 4°8 4°3 
tog [N(C-)/N(H-)] 529 S71 410 495 492 4°84 
Comp. Ilib spectral type M4(C7) M2(C6) K2(C3) Gs5(Co) Go 
sl log P 3°2 3°3 3°7 41 4°1 
— log [N(C~)/N(H>)] 343 389 224 2-32 2:32 








No. 3, 1958 Atomic and molecular negative ions in stellar atmospheres 263 


While carbon stars exist in which the C : H ratio approaches or even exceeds 
unity, there is no evidence that there are stars equivalently rich in oxygen. 
Further, the absolute photodetachment cross-section for O~- has been measured 
(Branscomb, Burch, Smith and Geltman 1958) and at 6000A is about 1/6 of 
the H~ cross-section. Accordingly, Table IV gives the relative O- concentrations 
and volume absorption coefficients [o(O-).N(O-)] at 6000A for the most 
favourable case of Ib (oxygen-rich) stars. Even with these extremely favourable 
assumptions, it is clear that O~ remains below the limit of detectability. 











TaBLe IV 
N(O-)/N(H-) in composition Ib dwarfs and giants 
Dwarfs Giants 

0) 200 1°70 1°40 1:20. 200 1°70 1°40 1:20 
Spectral type Ms M3 Ks-Mo Ka2 M4 M2  Ka2 Gs 
log P 50 5°0 50 4°95 3°2 3°3 38 471 
log [N(O-)/N(H-)] 2°81 2°65 2°65 2°62 2°96 2:88 2°75 2°66 
atl ot. ll 07009 0°006 0:006 0:006 O°013 O'OII 0°008 0006 
o(H).N (>) ’ . — 


* Taken at A6o00 A. 


Of the molecular ions of greatest interest, the OH~ ions will most likely 
affect late stars of normal composition, while CN~ and C,~ will be formed 
primarily in carbon stars. Table V gives the ratio of concentrations of OH- 
to H~ in main-sequence and giant stars of composition Ib. It is also instructive 
to compare the ratio of densities of OH~ and representative neutral molecules 
with extensive absorption bands in the visible and near infra-red spectrum, for 
example, CN. These data are also given in Table V. The OH~ concentration 
may be quite appreciable in Ms and Ss stars, in comparison with both the CN 
and the H~- concentrations, assuming the abundance of oxygen to be high 
enough. The absolute detachment cross-section for OH~ at A5500A is about 
one third that for H~ at the same wave-length. 

The ratio of the concentrations of OH~- to OH can be calculated by Saha’s 
equation alone, without knowledge of the composition, with the result that for 
dwarfs of any composition and of temperature T' = 2520 °K, log [pOH~)/p(OH)] 
=6-12. Thus molecular negative ions (like atomic ions) are only significant when 
the parent molecule has no continuous absorption of appreciable intensity in the 
observable wave-length range. 

The results for late-type giants of the same (Ib) composition (‘Table V) give 
somewhat smaller OH~ concentrations because of the reduced electron pressure 
and lower density. 


TABLE V 
N(OH-) in composition Ib dwarfs and giants 





Dwarfs Giants 
2) 200 1°70 61°40 1-20 2°00 1°70 =—-°40 1-20 
Spectral type Ms M3 Ks-Mo K2 M4 M2 Kz2_ G5 
log P 5O0 8650 550 4°95 32 33 38 41 
log [N(OH~)/N(H-)] 1°13 2°90 2°41 3°51 116 2°50 93°53 4°73 


log [N(OH-)/N(CN)] T'70 2°53. 3°10 5527 2,990 4°88 5:17 6:19 
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For carbon stars we first calculate in Table VI the concentrations of CN- 
in comparison with H~, CN and C, in cool stars of composition IIIb. In the 
calculation of p(C,) the values given by de Jager and Neven (1957) are used 
for p(C,°I1,,), and are affected by the error in the statistical weights. 





TaBLe VI 
N(CN-) in composition IIIb dwarfs and giants 
Dwarfs Giants 
(0) 2°00 1°70 1°40 1-20 2°00 1°70 3=—:«1*40—Ss«éd'*20 
Spectral type Ms M3 Ks-Mo Kz Ms5s(C7) M2(C6) K2(C3) Gs(Co) 
log P SO 50 50 4°95 32 8 63°3 38 41 
log [N(CN-)/N(H-)]) 3°34 1°99 O80 1°97 2°54 1°47 047 = 1°63 


alog [N(CN-)/N(CN)]} 3708 418 5:32 6°85 538 658 652 6:55 
Plog [N(CN-)/N(Ca)az] 3°02, #15 5°42 532 532 Hr 885-35 
a, Independent of composition 


b, Using de Jager and Neven’s figures for p(C.), which assume “II ground state and a 
statistical weight factor for C, of 9 (which is too large by the factor 3). 


We see from Table VI that CN- dominates H~ in all cases later than K2 
dwarfs (though no dwarfs of this composition are known at present) and G5 
(or Co) giants. 

An additional comment should be made about the most abundant molecule 
in late stars, H,. Wildt (1957) has discussed the possibility that H,~ ions 
contribute to the continuous absorption in low-temperature atmospheres. 
However, H, has a very large negative electron affinity (—3-58eV) and H,~ 
either has a very shallow potential well, or has a repulsive ground state and 
does not exist at all. The most detailed calculations of the H,~ potential function 
are those of Dalgarno and McDowell (1956), and they find EA(H,)= 
—3°58eV, and D(H,~)=o-15eV. Thus it is doubtful whether H,~ exists at 
all, and in any case it could not exist in stars in thermal equilibrium with Hg. 
Similar considerations apply to He~, which has a metastable, bound excited 
4P.. state, but of large positive total energy. 

3. Source of violet absorption in carbon stars.—It is of particular interest to 
investigate the possibility that CN~- contributes to the strong continuous 
absorption in the violet spectra of late N stars referred to in the introduction. As 
will be discussed below, this absorption is usually attributed to C, Swan bands 
(*II,,—*II,), to continuous absorption by C; and to the 4050 A Swings bands of 
C;, though SiC, may also play a part (Herzberg 1957). For this reason we 
would like to compare the absorption coefficient of CN- with both the 
unresolved Swan bands and the C, continuum. 

Let us first consider the Swan bands. In the calculations of de Jager and 
Neven, X*II,, was taken to be the ground state of C,, and only the A®II and B*Il 
states were considered in addition to the ground state in calculating the 
dissociation constants for C,. The possibility that the ground state of C, was 
the a!Z, state, and not X*I],,, was raised by Norrish, Porter and Thrush (1952). 
However, subsequent work by these authors (Norrish, Porter and Thrush 1953) 
suggested that the triplet state was the ground state, but that the singlet state 
lay less than a volt above it. And although Phillips (1957a) reported preliminary 
work by Brewer suggesting that *[], lay 1-2 eV above a!X, a recent investigation 
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of this point in Brewer’s laboratory by Hicks (1957) shows that X*II,, is the lower 
of the two states, with al having an excitation energy of only 0-34+0-17eV. 
De Jager and Neven were therefore justified in their choice of ground state, and 
their neglect of the a!X state in the calculations introduces very little error 
because of the low weight (0-5) of the state. Calculations on this basis of 
N(CN-)/N(C,3I1) are given in Table VI. The proper identification of the 
position of the X*J],, state of C, is important, because if it were true that X*I1, 
is excited by 1-2eV, the intensity of the Swan bands in absorption would be 
lower by a factor of about 250 (at @ =2) than would be concluded from the data 
in Table VI. 

Hicks (1957) also gives the absolute oscillator strength of the (0,0) Swan 
transition as 0-013. Using the relative transition probabilities given by Phillips 
(1957b) we can estimate the continuous absorption cross-section equivalent to 
unresolved bands of the Av=2 sequence at A4382A. ‘Thus the oscillator 
strengths for (2,0), (3,1) and (4,2) are about 0-39x10-%, 1-3x 107% and 
2-6 x 10-* respectively. Multiplying these f values by the ratio of the Boltzmann 
factor for the lower state to the vibrational partition function, and summing over 
the transitions in the sequence, we find 


[> @) 2 co N fo.” 
> foc’, v”+2)dA= — (10-8?) > f(v’,v” +2) 2 e ) ~2x to-*omtA. 
0 vw” =0 





” 


v= 


4 


The sequence extends over a spectral range of about 100 A, so that the equivalent 
continuous absorption cross-section is about 2x 107%cm?. The continuous 
absorption coefficient of CN~ is unknown, but may be 10 or 100 times greater. 
Assuming the factor to be 50, the concentration ratios of Table VI correspond 
to ratios of absorption coefficients of CN- to C,(Av=2) of 10~% for C7 giants 
and 5 x 10~-* for Ms dwarfs of the assumed composition. Thus we expect the 
Swan band absorption at 14380 to exceed the CN~ continuum by a factor of 
10? — 10%, subject to the assumptions made. 

The rapidly decreasing oscillator strengths with increasing Av in the Swan 
bands suggests that the Av=3 sequence near the peak of the continuous 
absorption in N stars will be weaker than Av=2 by a factor of about 5 at O=2. 

Perhaps a more significant comparison is between the CN~ continuum and 
the C, continuum. Calculation of the C, concentration is handicapped both by 
its sensitivity to the atomic carbon concentration and by the unavailability of 
accurate free energy functions for C;. However, Goldfinger (1955) has 
estimated the free energy function (F°—£,°)/T for C, to be about 66-5 +3 at 
2400 °K from the properties of similar tri-atomic molecules. If we use this 
figure and compute dissociation constants for C3, we find that p(C;)/p(C,) is 
small for @ < 1-4 and is of the order of 1 for © = 1-7 in giants of composition IIIb. 
To calculate p(C,) at ©=2 would require a recalculation of the chemical 
equilibrium given by de Jager and Neven. An adequate approximation for the 
late carbon stars may be found in the assumption that p(C;) is equal to the value 
of p(C,) given by de Jager and Neven. For © < 1-4 this will be an over-estimate. 

Phillips (1957a) states that the excitation potential of the lower state of the C, 
continuous absorption is at least 2-6 eV. In 'Table VII we combine this minimum 
result with the assumption that p(C,)~p(C,;) to find log [p(CN~-)/p(C,*)] 
where the asterisk refers to the excited state of C, from which the continuous 
absorption occurs, Even if we take one-third of the total carbon content of 


18 
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the atmosphere as the value of p(C;), the ratios given in the table are reduced 
by only a factor of 30 for@=2. ‘Thus, although this calculation is quite crude, 
we may expect CN~ to be a significant competitor to C; continuous absorption, 
if Phillips is correct in his assignment of a 2-6 eV excitation energy to the lower 
state of the C; transition. As Phillips himself pointed out, however, the optical 
data on which the figure of 2-6 eV is based are in conflict with the known heats 
of sublimation of C, C, and C, (Brewer 1958), and recent work, privately 
communicated by Brewer, strongly suggests that the laboratory emission of C, 
observed by Phillips and Brewer (1955) is not a true continuum but an 
unresolved complex band system of low excitation potential. The absence of 
this emission from comets (in which the Swings bands are strong) suggests 
that the excitation potential is not zero, though the new laboratory data make it 
unlikely that it is as high as, say, 2 eV. Hence it is somewhat doubtful 
whether CN- is sufficiently abundant in stellar atmospheres, compared with 
excited C,, to contribute appreciably to the violet absorption in N stars. 


Tasie VII 


Comparison of CN- absorption in IIIb dwarfs and giants with the C3 continuous absorption 
using Phillips’ data on the states of C3. 





Dwarfs Giants 
1) 2°00 1°70 1°40 1:20. 200 1°70 1°40 1:20 
Spectral type Ms M3 Ks-Mo Kz Ms(C7) M2(C6) K2(C3) Gs(Co) 
log P 5O 86570 50 4°95 32 3°3 3°8 4'1 
alog N(CN-)/N(C;*) 2°22 oO°57 1°06 2°44 0°52 1°03 2°52 2°47 


a, Assuming that the lower state of the C, continuum is excited by 2:6 eV, in accordance 
with the work of Phillips and Brewer (Phillips 19574). 


If the lower state of the true C, continuum is 5 eV above the ground state of C,, 
then the calculated CN~ concentration in late carbon stars will exceed the C,* 
concentration by an enormous factor, but this high excitation would probably 
also prove that the C; system observed in the laboratory cannot be the continuum. 
We do not suggest of course that the laboratory spectra in King furnaces and 
acetylene-oxygen flames should be associated with negative ions, for the electron 
pressure in these sources of the continuous emission (Phillips and Brewer 1955 ; 
Kiess and Broida 1956; Marr and Nicholls 1955) is far too low. One might 
be forced to find another explanation for the laboratory ‘‘ continuum ”’ and it 
would then remain to be proven that the CN~ in stars was competitive with the 
new mechanism. The observed similarity between the laboratory continuum 
and the continuous absorption in carbon stars (McKellar and Richardson 1955) 
is still strongly suggestive of a common origin. Perhaps one might interpret 
our calculations of the CN~/C,* ratio in stars as evidence against a true C, 
continuum, rather than as specifically supporting CN~- as the absorption 
mechanism in the stars, particularly in view of the observational evidence 
mentioned below. 

It is interesting to note that the wave-length dependence of the absorption 
curve deduced by McKellar and Richardson (1955) for Y CVen, for example, 
closely resembles not only the C, curve observed in emission by Phillips and 
Brewer, but also the curve which one would expect for CN~- absorption with 
EA(CN)=2:92 eV, The main maximum at 4000 A and the subsidiary maximum 
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at 4300 A would then correspond to the peak absorptions for the Av=o and 
Av=—1 sequences, respectively. The apparent thresholds of the two maxima 
imply w,(CN~)~ 1000 cm~, which is of the correct order of magnitude, although 
it is about half the vibrational energy for the ground state of CN. However, 
the correlation between the continuous absorption and the strength of the CN 
bands is generally poor; thus, since the ratio p(CN~-)/p(CN) varies by a factor 
of only 6 throughout the temperature range assumed for giants, considerable 
variations in electron pressure must be present if our identification of CN- is 
correct. 

In this connection, one should also note the observation by Feast (1955, 1957) 
that the strong continuous ultra-violet absorption is present in the peculiar 
stars AM Cen and GP Ori, which do not show C, bands and are believed to be 
intermediate in C: O ratio between N and M stars. In view of the great strength 
of the sodium D lines in both these stars and the N stars, Feast suggested that 
the ultra-violet opacity might be due to unresolved NaH bands. On the other 
hand, it is possible that high Na abundance might increase the electron pressure 
above the value we have assumed, thereby favouring the formation of CN-. 
In any case, one must reckon with the possibility that more than one species 
contributes to the ultra-violet continuous opacity in N stars. 

The concentration of C,~ can be expected to be rather similar to that of 
CN-, using the IIIb composition of de Jager and Neven, for the electron affinity 
of C,, while poorly known, is similar to that of CN (Honig 1954). However, 
the expected E** threshold energy dependence of the detachment cross-section 
will make C,~ a less important ion in the wave-lengths in which observations 
can be made. 


4. Further relevant features in stellar spectra. 


4.1. Sun and “‘ordinary’’ late-type stars.—The spectral energy distribution 
and limb darkening in the ultra-violet continuous spectrum of the Sun has 
been studied by Canavaggia and Chalonge (1946, 1950) down to 3150 A. From 
these studies, Chalonge, Divan and Kourganoff (1950) concluded that there 
must be a strong additional source of absorption in the ultra-violet, apart from 
H and H-. A detailed study of the properties of this unknown absorption 
‘*X”’ was carried out by Michard (1953), who concluded that the additional 
mass-absorption coefficient increased with diminishing wave-length from a 
threshold in the neighbourhood of As5ooo A, and that its value at 43200 A increased 
with depth to a maximum of about 0-37 at optical depth 7,,,.~0°55 and then 
diminished with increasing depth. Michard found that X could not be accounted 
for satisfactorily by metal continua and pointed out that its behaviour was con- 
sistent with Swings’ suggestion that negatively ionized molecules play a part. 

Two features of X seem to lend confirmation to Michard’s conclusion. The 
first is the behaviour with wave-length, which turns out to be very close to an 
E32 law. An inspection of Table I immediately suggests C,~ as a possibility. 
The second is the behaviour with depth; the continuous absorption of an atom or 
negative atomic ion with an electron affinity of only 2-5 eV (corresponding to 
a threshold at A5000 A) would increase monotonically with depth in any reasonable 
model solar atmosphere, as would the effect of an underestimate of the cross- 
section of H-. To obtain the decrease with depth for 7; 9)>0°55, one needs 
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a molecule with a dissociation potential of at least 3 eV, neglecting competition 
effects. 

Michard’s analysis of the additional absorption is open to objection because 
of the neglect of the metallic continua. His data for relative absorption co- 
efficients at different wave-lengths were therefore re-analysed for the layer with 
T5000 = 0°40, 7’ =5750°K, assuming the absorption at A5oo0 A to be given entirely 
by known sources as included in Vitense’s tables (Vitense 1951). The two- 
thirds power of the residual mass absorption coefficient, «,..—*;caio. WAS 
plotted againt 1/A and found to fit tolerably well to a straight line cutting the 
1/A —axis at A=4100+ 200 A (EA(C,)=3:0+0°2 eV). Unfortunately the first 
five points out of eight are all somewhat uncertain, the first three because the 
residual absorption is so small, and the next two because they are in the neighbour- 
hood of the Balmer limit, so that not too much significance is to be attached to 
this result in view of the probable effects of line absorption. Furthermore, 
the inclusion of the metallic continua led to a reduction by a factor of 4 in Michard’s 
values of k,.4;—*acacs although the variations with wave-length and with 
depth were nearly unchanged. 

Using Michard’s model of the solar atmosphere with de Jager and Neven’s 
computations of the C, abundance for composition Ia, assuming a dissociation 
energy of 6-7 eV for C,, it was found necessary to assume an impossibly large 
absorption cross-section of 2:5 x 10-14 cm? for C,~ at A3200 A to account for X 
in this way. As there does not appear to be any molecule with the appropriate 
properties and with sufficient solar abundance to account for the additional 
absorption in the solar ultra-violet, we must conclude provisionally that X is a 
combination of metal continua and line absorption. 

In normal stars of spectral type later than F5, colour temperature observa- 
tions by Woolley, Gascoigne and de Vaucouleurs (1954) suggest an ultra-violet 
depression of the same general type as is found in the Sun, increasing towards 
lower temperatures and higher luminosity. This result is in qualitative agree- 
ment with that of Keenan (1954) already quoted. In view of the luminosity 
effect and the results of our analysis of the solar data, we cannot conclude that 
this depression is due to negative ions. 

4.2. R Coronae Borealis.—The photographic region of this peculiar supergiant 
Fs5-8 variable has been studied at maximum light by Berman (1935) and at 
minimum light by Herbig (1949). Both the absorption spectrum and the 
emission spectrum that appears at light minimum are virtually free from hydrogen 
lines, while carbon features including those of C, are prominent in absorption. 
At minimum light Herbig found a broad emission feature degraded shortward 
from A3882, which he has tentatively attributed to the well-known Av =o sequence 
of CN. From his curve-of-growth analysis of the 4A4300-4900 A region at 
maximum light, Berman deduced an atomic composition with 69 per cent carbon 
and 27 per cent hydrogen. It follows (‘Table III) that C~ will be more than five 
times as abundant as H~ in its atmosphere, and that C~ is the dominant source of 
continuous absorption unless its photodetachment cross-section is very small. 
When laboratory measurements of this cross-section become available, it will 
be of great interest to construct a model atmosphere for R Coronae Borealis on 
this basis. It is also possible that careful spectroscopic study of this star will 
reveal the absorption threshold of the C~ ion and hence provide a more precise 
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figure for the electron affinity of atomic carbon than is possible in laboratory 
studies. 

4.3. S stars.—In addition to the well-known bands of TiO, CN, LaO, YO 
and ZrO, the visual and near infra-red regions of S stars contain numerous 
unidentified bands with heads at 45849 A and in the range AA8250-8820 A (Keenan, 
1957). ‘These bands are all degraded longward and thus can certainly not 
represent ionization or detachment continua. However, in view of the results 
obtained earlier for the OH~ ion, with a wave-length threshold in the neighbour- 
hood of A7000, it is of interest to examine the spectra of S stars in this region, 
since those of M stars (which probably have about the same oxygen abundance 
as S stars) are obscured by the TiO band at 7054 A. The TiO bands vary in 
strength among S stars, and there is thus some hope of finding suitable individuals 
in which the OH~ absorption (if present) might be detected. Spectra of this 
region with a dispersion of the order of 2000 A/mm have been obtained with the 
Case Schmidt telescope (Nassau and van Albada, 1948). Unfortunately the 
region still has numerous stellar bands of CN and atmospheric bands of H,O. 
The faintness of S stars makes it very difficult to carry out the detailed spectro- 
photometric study that will be required to detect the rather faint hypothetical 
dip in the continuum due to OH-. 

5. Conclusions.—We have shown that CN-, and perhaps C,~, may play an 
important part in the ultra-violet continuous absorption of carbon stars. On 
the other hand, our attempts to identify other negative ions have led to the 
conclusion that special circumstances must prevail for C~ and OH~ to be 
important. O- and O,~ are unlikely to be of interest in any celestial objects. 

The following types of investigation are needed to shed further light on the 
problems that have been raised: 

1. Accurate heterochromatic spectrophotometry of the visible and infra-red 
spectra of appropriate stars, particularly carbon stars, late S stars, and G and K 
giants deficient in hydrogen. High dispersion is needed to separate line and 
band absorption from continua, but the effect sought varies slowly with wave- 
length, which will make such studies difficult. 

2. More accurate determinations of the composition, total gas pressure and 
electron pressure of the specific objects under investigation. These are needed 
for the analysis of the observations. 

3. Laboratory measurements of absolute oscillator strengths and vibrational 
transition probabilities of the bands and continua of neutral molecules, to permit 
calculation of the competition between neutral and negatively ionized molecules. 
If the initial states of these transitions are excited, their statistical weights and 
excitation potentials must be determined. 

4. Actual observations of the absolute cross-sections for photodetachment 
for those astrophysically interesting ions which have not yet been studied in the 
laboratory. Sucha programme is under way at the National Bureau of Standards. 

We are grateful to L. Brewer, M. W. Feast, G. Herzberg and A. McKellar 
for their helpful criticism. 


University College, Royal Greenwich Observatory, 
London: Herstmonceux Castle. 
1958 March. 
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ON THE REDUCTION OF POLARIZATION MEASURES 
AFFECTED BY INSTRUMENTAL POLARIZATION 


A. F. Wesselink 
(Communicated by the Radcliffe Observer) 


(Received 1958 March 10) 


Summary 


The effect of oblique reflection on measures of polarization has been 
studied. Simple formulae are derived whereby observations of stars whose 
polarizations are known can be used to calibrate observations subject to 
oblique reflection. The method is illustrated by Mrs E. van P. Smith’s 
observations. The utility of the quantities pcos 20, psin2@ for the 
discussion of internal accuracy and comparison of independent series is 
stressed. 





It is well-known (1) that in the measurement of the polarization of stars due 
attention must be paid to various possible sources of systematic error. So far, 
two effects have been noted by observers. First, it has been found (1) that the 
sensitivity of the inclined photocathode of the RCA multiplier varies with the 
position angle of the linearly polarized light. Even the head-on cathode of the 
EMI tube shows a slight but significant position angle effect. This difficulty, 
however, can be overcome conveniently by observing stars for which the 
polarization can be expected to be absent or negligible. 

Another instance of instrumental polarization occurs when the stellar light 
passes the analyser after oblique reflection. If linearly polarized light is reflected 
obliquely from a metal surface, the component of the vibration in the plane 
containing the normal is reduced differently from the component at right angles 
to this. 

Although coherence remains preserved, the two components suffer unequal 
shifts in phase. As a consequence, the linearly polarized part of the stellar light 
is transformed into elliptically polarized light before passing the analyser of the 
photometer. The subsequent complication in deriving the true polarization 
from measures thus distorted induced observers to avoid carefully oblique 
reflection when measuring polarization. 

The only observer, to the writer’s knowledge, who successfully observed 
stellar polarization, though the optical column contained oblique reflection, has 
been Mrs E. van P. Smith (2). She used the 60-inch reflector of the Boyden 
Station at Bloemfontein in her noteworthy study of polarization in southern B 
stars. The primary mirror of this telescope does not contain a central hole and 
the telescope can therefore be used only in the Newtonian focus. The Newtonian 
flat is the oblique reflecting element in the optical column. Many stars were 
measured with the Boyden 13-inch telescope as well; since the results with this 
telescope may be considered reasonably free from systematic error she was able 
to reduce the 60-inch observations to true polarizations. In these reductions 
Mrs Smith used graphical methods. 
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It is the purpose of the present article to show that, as an alternative to 
Mrs Smith’s procedure, simple formulae can be derived whereby observations 
of stars whose polarizations are known can be used to calibrate observations 
subject to oblique reflection effects. 

If p is the degree of polarization, 6 the position angle of the electrical vector, 
and « the position angle of the analyser, we have for the intensity of the light 
as measured with a polarimeter devoid of instrumental polarization for the 
intensity as function of «: 


i(a) =1+pcos2(0—«) (1) 
=1+ pcos 20 cos 2«+ psin 20 sin 2«. (1') 
If we put: 
x=p cos 20 (2) 
y=psin20 (3) 
we have 
i(x)=1+%cos2a+ysin 2a. (1”) 


In an unpublished theoretical investigation, I have found that the intensity 
of the light when recorded after oblique reflection can be similarly expressed: 


i’(a)=1+.' cos2a+y’ sin 2a. (4) 


In (4) the quantities x’ and y’ are each simple functions of x and y that contain 
constants depending on the condition of the mirror, the angle of incidence, and 
the wave-length of the light. The theoretical relations just mentioned are very 
nearly linear between the quantities x, y, x’ and y’ for the small polarizations 
actually occurring among the stars. 

The conditions in a telescope are less simple than those assumed in the theory : 
for instance there is a considerable spread in the angles of incidence of the various 
rays constituting the converging beam. Though there are undoubtedly assump- 
tions in the theory that are realized only approximately in practice we may hope 
that linear reduction formulae of the following form represent the observations 
satisfactorily : 

x=ax' +by'+c (5) 


y=dx' +ey' +f (6) 


in which the constants have to be found from the measures x’ and y’ of stars 
with known polarization x and y. If we call p’ and 6’ respectively the distorted 
polarization and position angle and define them similarly to x’ and y’ as p and @ 
are related to x and y by the formulae (2) and (3), we see that a possible error in 
the assumed zero point of the analyser will be absorbed automatically in the 
constants of the formulae (5) and (6). 

In order to show the applicability of the formulae in actual practice we applied 
formulae (5) and (6) to the observations of Mrs Smith who fortunately gives the 
60-inch results p’ and @’ and the 13-inch results p and 6 for a number of stars. 
The Bloemfontein 60-inch polarimeter was provided with a flat mirror whose 
function was to reduce the oblique polarization caused by the Newtonian flat. 
This circumstance of double oblique reflection certainly affects the numerical 
values of the constants in the formulae (5) and (6) but should not necessarily 
influence the form of these relations. 


No. 3, 1958 measures affected by instrumental polarization 273 


Mrs Smith divides her 60-inch results in two groups, (a) and (b), each of 
which may be considered homogeneous with respect to instrumental polarization. 
For each group, two least squares solutions were made of the forms (5) and (6), 
the results of which are: 


( x= +0°899x' + 0-412y’ — 0°0275 
+0:022 +0°026 +0°00II (7) 
(a) 33 stars { 
y= —0°375x' +0°871y’ + 0°:0220 (8) 
{ +0015 +0°018 +0°0008 
[x= +0:904x' + 0:380y’ — 0°0428 
+0°015 +0°01I9 +0°0008 (9) 
(b) 39 stars 3 
y= —0°413x’ + 0-g20y’ — 0:0028 (10) 
+O°OII +0°014 +0°0005. 





It should be remembered that the polarization in these formulae is expressed 
in the manner nowadays customary in astronomy, i.e. p is the difference in 
magnitude between the light at «=6@ and at «=6+g90°. For small polarizations 
the percentage polarization is 46-05 times larger. If the formulae (7)-(10) are 
rewritten with percentage polarizations there will be no change in the coefficients 
of x, y, x’, y’ but the constant terms have to be multiplied by 46-05. 

The success of the reduction method in this particular example can obviously 
only be judged from the absence of a systematic run of the residuals with various 
parameters. I have therefore plotted these residuals in the solutions for x against 
x’, y’, and x, those of the solutions for y against x’, y’ and y. 


TABLE I 


Number of sign alterations Q in residuals of solutions when arranged according to vartous 
parameters. Qe is the expected number on the assumption of complete randomness 


Parameter Q Qe Residual in 
x’ 12 (16) 
y’ 15 (16) \ x, equation (7) 
x 16 (16) 
i 15 (16) 
y’ 17 (16) k y, equation (8) 
y 18 (16) 
x” 17 (19) 
y’ 15 (19) } x, equation (9) 
x 19 (19) 
x’ 19 (19) | 
y’ 24 (19) j y, equation (10) 
y 16 (19) 


Table I gives the number of changes in sign in each plot. The expected mean 
number of such changes is }(m—1) when n is the number of observations or 
16 for the (a) solutions and 1g for the (4) solutions. It is seen that this sign test 
is adequately satisfied in all cases. 
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Residual from equation (7) 


Residual from equation (8) 


The plots are graphically shown in Fig. 1 only for the (a) observations. 
The formulae (7)-(10) can also be written as follows: 


x= +0:989p’ cos (26’ — 24°*7) —0°0275 (7 
y= +0°949p’ sin (20’ — 23°-3) +.0°0220 (8 
x= +0°982p’ cos (26’ — 22°-8) — 0°0428 (9 
y= +1-008p’ sin (20’ —24°-2) — 0-0028. (10 


‘) 
‘) 
‘) 
‘) 


It is seen that the constant phase angle is sensibly the same in all four formulae 
(24°). The simplest explanation for this is that there might have been an error 
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Fic. 1.—The figure illustrates for the (a) observations the behaviour of the residuals of the 
solutions with respect to various parameters. The top three diagrams give the residuals of equation 
(7) (ordinates) as functions of respectively x’, y’ and x whereas the bottom three diagrams exhibit 
the residuals from solution (8) as functions of x’, y’ and y. The scales in the abscissae and the 
ordinates are equal. Each diagram is 0™-06 high and o™-30 wide. The equation is indicated 
at left whereas the parameter is inserted at the right. 


x’ 
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of 4x 24°=12° in the zero point of the analyser, though other explanations 
cannot be ruled out entirely. 

We have seen the utility of the quantities x and y for the reduction of 
polarization measures on to a standard system. They are also more useful than 
p and @ for the study of internal accuracy or whenever different independent 
series of measurements are to be compared. 

For, whenever polarization is measured by measuring brightness with an 
analyser set subsequently in a number of position angles (m) evenly distributed 
between 0° and 180°, it can be shown that the uncertainties in x and y are exactly 
the same and equal to the uncertainty in p which is +eV8/n where « is the 
uncertainty in the determination of a single magnitude. The useful feature is 
that the uncertainty is only dependent on the accuracy « and the number of 
measures, but is independent of the value of the polarization itself; by contrast 
the error in the position angle is «/pV 2/n which for small polarizations may 
become large. Though the Bloemfontein measures, discussed in this article, 
were made in a slightly different manner, we do not expect large differences in 
weight to have resulted in x and y. Therefore the least squares solutions were 
made giving equal weight to the equations of condition throughout. From the 
solutions I find a mean error in either x or y of +0™-005. Hence the mean error 
in p is also +0™-o05 independent of the amount of polarization whereas the 
mean error in the position angle is correspondingly computed as + 0:0025/p radian 
or +o°15/p degrees. 


Radcliffe Observatory, 
Pretoria : 
1958 March 3. 


Appendix.—Though not connected with the problem of eliminating the 
effects of oblique reflection effects, the general usefulness of the quantities x and y 
has been discussed before by Tripp (3), whose results are in agreement with those 
given in the present paper. A theoretical study of the effects due to oblique 
reflection by metal surfaces was made by Seares (4) in a paper entitled ‘‘ The 
displacement curve of the Sun’s general magnetic field’’. 

I am indebted to a referee for drawing my attention to both these articles. 
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Summary 

A method is described whereby the amplitude and phase of the complex 
Fourier transform of a spatial brightness distribution of small angular extent 
may be uniquely determined from a series of measurements with a triple 
interferometer system. Absolute measurement of the amplitude function is 
available, whilst measurements of phase are relative to a datum obtained 
at short aerial spacings. A practical radio frequency interferometer incor- 
porating the principle is described and its operation is discussed. 





1. Introduction.—Unique measurements of the brightness distribution across 
the cosmic radio sources have hitherto been limited to those sources which 
subtend a large angular diameter. It has been possible to obtain orders of 
magnitude for the angular diameter of some of the smaller intense sources and 
in a few cases to give models for their general intensity distributions (Hanbury 
Brown, Jennison and Das Gupta 1952, Mills 1952, Jennison and Das Gupta 
1953, 1956). ‘These models have been obtained from the amplitude functions 
of the Fourier transforms of the source brightness distributions without corre- 
sponding information with regard to the phase. As such, the solutions are not 
unique and may not correspond in detail to the actual fine structure of the sources. 

The amplitude function of the Fourier transform of the brightness distribution 
across a source may be determined from the visibility of the fringes in the reception 
pattern of an interferometer. The phase of the transform may be obtained from 
the phase of the fringes relative to a known datum. This datum may be obtained 
by accurately measuring the position of the source over various interferometer 
baselines. 

The practical limitations imposed on these measurements are severe and 
render conventional techniques untenable when the baseline of the interferometer 
system exceeds about one thousand wave-lengths. At greater distances, when 
it is necessary to use a radio link between the channels, the measurement of 
fringe amplitude may be accomplished with an uncertain error if a powerful 
automatic gain control system is used on the receiving equipment to maintain 
the mean level of the signal constant over the period of the observations. The 
aerial gain, noise factor, bandwidth and dispersion are assumed constant or 
must be continuously calibrated (Mills 1953). If both the receiving aerials are 
sufficiently large the uncertainty may be removed by dispensing with the automatic 
gain control and measuring the increment in total power due to the source in each 
channel (Jennison and Das Gupta 1956). Neither of these methods allows for 
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error due to phase dispersion in the systems and the latter method has the practical 
disadvantage of requiring a very large portable aerial. The method used by Smith 
(1952) is most reliable on short baselines but is not applicable where the distance 
between aerials is considerable. In these circumstances intermediate checks 
cannot rapidly be made and the gain of the most distant station cannot necessarily 
be relied upon to remain constant over long periods or on removal to a new site. 

The conventional method for the measurement of the phase of the fringe 
system by timing the transit of the central fringe is not applicable at long baselines, 
especially when radio links are used to span the baseline. The combinations 
of variable phase errors in the electronic equipment and in the local oscillator 
and intermediate frequency signal link paths, together with the difficulties 
arising from unequal signal paths when the terrain is not level, render the phase 
of the transform indeterminate at spacings in excess of about one thousand 
wave-lengths. 

An interferometer system designed to overcome the above defects was 
constructed in 1954 and its principle of operation is discussed below. 

2. Principle of operation.—The main features of the apparatus are shown in 
Fig. 1. It consisted of three aerial systems, tuned to a frequency of 127 Mc/s, 
spaced in a line along which the measurements were to be made. One of the 
aerials was a large fixed array whilst the other two were small and portable. 
Associated with the portable systems were radio link transmitters to relay the re- 
ceived signals back to the fixed station. Provision was made to rotate, continuously 
and synchronously, the phase of the two relayed signals prior to equalizing 
the time delays, and combining all three signals in pairs in three switched inter- 
ferometers, the outputs of which were displayed on a three-pen recorder. The 
total power received from the fixed aerial was displayed on a separate recorder. 


The measurement of phase 


The relative phases of the signals induced in the three aerials are indeterminate 
by the time they reach the switched interferometer units, due to the combinations 
of errors previously referred to. Thus the fringe systems produced at the outputs 
of these units and displayed upon the pen recorder are also of uncertain phase 
relative to an absolute datum based upon the time of transit of the source. It 
will be shown that the relative phase contribution due only to the brightness 
distribution of the source may be obtained by a comparison of the arguments 
of these fringe systems. 

Treating the signals from aerial A as the phase reference and taking the 
origin of the transform at zero within the source, let &,,, gc and €,, represent 
the phase of the transform at spacings AB, BC and AC respectively; let w,(t) 
and w,(t) represent the time variable phase rotation of the scan at B and C 
respectively; let ys, and ys represent the phase angle introduced by the position 
of the source in relation to the collimation plane at B and C respectively; and 
let 5, and 5, represent the phase error in the equipment at B and C respectively. 

The argument of the fringes when channel B interferes with A will be 


Ean t+4,+8,+o,(t) (1) 
and the argument of the fringe when channel C interferes with channel B 
= fot (po—Y1) + (82 —8,) + w9(t) — w, (2). (2) 


The argument of the fringes formed by channel C interfering with channel A 
= Exot pet. +w,(t). (3) 
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Adding the arguments of the two fringe systems AB and BC, we obtain 


Eat pot pe +82 +w,(t). (4) 
Comparing equations (3) and (4) it will be seen that they represent the arguments 
of two fringe systems of identical frequency and relative phase 


$=£xo— (Exam +c): (5) 


This result is independent of the phase errors and the angular velocity of the 
phase rotation. It enables an estimate to be made of relative phase contributions 
due to the source structure only, and is applicable in the general case where the 
aerials are spaced by arbitrary distances. It will be observed that the position 
of the source relative to the collimation plane does not appear in the expression 
Exsc— (€4n+&gc) and hence absolute phase is not directly obtainable. 
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Fic. 1.—The principal features of the phase-sensitive interferometer. 
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It is apparent that if the spacings AB and BC both lie within the first maximum 
of the Fourier transform of a simple symmetrical distribution, then €,3 =£,,=0. 
In these circumstances $=€,,.=0 if AC lies within the first maximum and 
¢=£,c=7 if AC lies within the second maximum. 

If the transform is complex but the positions of the aerials are adjusted so 
that AB= BC, equation (5) reduces to 


$=fsc—2€ap 
and a direct check on the relative phase over twice the spacing may thus be 
obtained. 

In the more general case, when the transform is complex and the aerial 
spacings unequal, the expression (5) may be used for the measurements in the 
following manner. 

The two portable aerials B and C are moved away from the home station A 
with their baselines in any ratio until a phase displacement is recorded. This 
may be most simply done by keeping one station, B, at some convenient distance 
from A and moving the other until a displacement is observed. Within this 
range of spacings we may put 


€,n=pco= constant. 


It is apparent that the phase function does not vary, within observational error, 
over intermediate spacings within the maximum separation AC, and hence both 
portable stations may be moved to more distant sites B’ and C’ subject only to 
the limitation that AB’ < AC and B’C’<AC. 

The relative phase occurring over the extension of the baseline to points 
in the range AC’ may now be established. The process may be repeated until 
the whole phase function is eventually mapped. It will be observed that this 
general method gives rise to a cumulative error, and the accuracy of the 
determination of the transform depends upon the ratio of the maximum aerial 
spacings used to the spacings at which a displacement is first observable, and the 
error in the assessment of the phase at this spacing. 


The measurement of amplitude 


From direct measurement on the pen recordings it is possible to form the 
function 7 such that 
_ [4B] x AC] 
1 [BC|x [AP ’ 
where |AB|, |BC| and |AC| are the moduli of the fringe systems between 
channels A and B, B and C and A and C respectively, and |A|? is the total power 
due to the source recorded in channel A. 

It will be shown that the function 7 is a measure of the relative amplitude and, 
in a special case, the absolute amplitude, of the Fourier transform of the intensity 
distribution of the source. 

Let g,«, go8 and gsy be the mean amplitudes of the signal at the outputs of 
channels A, B and C, where g,, go, g3 and «, B, y refer to the overall voltage gains 
and input signal levels respectively. Let f(ab), f(bc) and f(ac) be the amplitudes 
of the transform for baselines AB, BC and AC. Then 


= |AB| x |AC| es G12,%g8f(ab) x Gag ,agsyf (ac) - f (ab) x f (ac) ‘i G,G, 
1 [BC[x [A] Gag aPeayf (be) x G4(g,x)? f (bc) CG. , m 
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where G,, G, and G;, are the conversion gains for the cross multiplying recorder 
systems between channel A and B, B and C, and A and C respectively and G, is 
the conversion gain of the total power recorder system in channel A. The 
expression for » in equation (6) is not dependent on the gain of any part of the 
equipment other than the final cross multiplying and total power recording 
units. The equipment was designed so that these units were operated at high 
level and were extremely stable. Sockets were provided so that, immediately 
prior to and following a run, the inputs normally fed from the separate channels A, 
B and C could be connected together and provided with a single noise signal. 
From the constant deflection of the recording pens this yielded a normalizing 
function : 
G,G, 

“2. (7) 
corresponding to a correlation coefficient of unity between the channels. 

Substituting the normalizing function K in equation (6), the general 
operational function y’ may be obtainable, where 


y= -ABIxIAC]_ _ flab) xflac) 
K x [BC] x JAP fe) 


K 





(8) 


This expression only involves terms in the amplitude function of the Fourier 
transform of the brightness distribution across the source. It will be observed 
that in the special case where station B is midway between stations A and C, 
f(ab)=f(bc) and equation (8) reduces to 


9’ =f (ac). 
The absolute value of the amplitude of the transform may therefore be determined. 

In general, when AB# BC, equation (8) may be used to obtain a relative 
measurement of the amplitude at one spacing from a knowledge of the amplitude 
at the other two. Measurements may therefore be made concurrently with the 
phase measurements, commencing at short spacings and extending the baselines 
as each value of the amplitude function is established. 

3. A practical radio frequency interferometer.—The principles outlined above 
were incorporated in an interferometer constructed at Jodrell Bank and operating 
on a frequency of 127 Mc/s. 

The aerial A consisted of 160 full-wave dipoles in a collinear broadside array. 
This aerial had a beam width of approximately 3°x10°. Aerial B usually 
consisted of 20 full-wave dipoles in a collinear broadside array giving a beam 
pattern 24°x10°. Aerial C was either a single folded dipole with reflecting 
screen or an array of four similar dipoles. 

The basic lay-out of the electronic equipment is shown in Fig. 2. The out- 
stations each employed two crystal oscillators to feed their first and second mixers. 
The beat between these two oscillators was transmitted as a coherent local oscillator 
reference. The noise signal was transmitted from the broadband class A amplifier 
following the second mixer. Transmitter powers were normally about 0°25 
and 1 watt from the signal and local oscillator respectively. The relay transmitting 
and receiving aerials consisted of single dipoles with parasitic reflectors, the 
polarization of the signal and local oscillator aerials was crossed and orthogonal 
to the corresponding aerials of the other radio link, 
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The receiving equipment for the radio links consisted of separate preamplifiers 
for the local oscillator and signal channels. The low level of the signal from the 
received local oscillator was converted coherently to a signal of the order of volts 
by being mixed with the tenth harmonic of a master local oscillator on 11-5 Mc/s, 
and subsequently passed through a narrow band amplifier at approximately 
30 Mc/s. The output of this unit was applied to a mixer following the signal 
link preamplifier and yielded an output on an intermediate frequency of 12 Mc/s. 





























* ee: 
i © 


Fic. 2.—Basic block diagram of the phase-sensitive interferometer. The numbers in the 
blocks correspond to: 1, amplifier; 2, buffer; 3, mixer; 4, oscillator; 5, frequency multiplier; 
6, delay line. 


The corresponding equipment in channel A consisted only of a preamplifier, 
a mixer combining the 115 Mc/s master local oscillator frequency with the 
signal frequency, and an I.F. amplifier on 12 Me/s. 

The 12 Mc/s I.F. amplifiers in all three channels were followed by mixer 
valves converting the frequency to 500kc/s. The conversion was performed 
by mixing with the master local oscillator on 11-5 Mc/s, directly in the case of 
channel A, but via continuously rotating phase shift networks in the case of 
channels B and C, These phase shift networks consisted of resistance capacity 


19 
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quadrature circuits applied to a four pole condenser system having a rotating 
search vane of high dielectric constant. The amplitude of the output of this 
unit was substantially constant, and the tracking error in the phase angle was 
better than 2°. The phase rotating systems in channels B and C were mechanically 
coupled to friction rollers adjustable in position on a single conical shaft driven 
from a synchronous electric motor. A large but constant ratio between the rates 
of phase rotation of the two channels was thus available. 

The three signals at 500 kc/s were confined to matched bandwidths of 200 kc/s, 
or, in certain measurements, 50 kc/s. The signals from channels A and B were 
then delayed by suitable lengths of artificial lines before being applied to the 
cross multiplying and total power recording unit. 
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Fic. 3.—Block diagram of one of the cross-multiplier units shown as a cross in Figs. 1 and 2. 


The three cross multiplying circuits and the total power recording circuits 
were identical except that the two input leads of the total power recorder were 
connected together. A basic block diagram of a single section is shown in Fig. 3. 
500-cycle square waves were used for switching, whilst the actual multiplication 
was performed in an octode valve; an auxiliary circuit with a balancing valve 
eliminated single channel noise by cancelling the effect of non-linearity on the 
variable-» signal grid. Two phase splitters and switching circuits only were 
required to serve the three multipliers and one square law detector which operated 
on a similar principle. 

4. The equipment in operation.—The equipment was extensively tested on 
short baselines within the range of readings previously obtained with conventional 
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techniques. During the tests a slight dispersion in the radio links was observed, 
This was reduced by improved design and finally removed entirely by the 
introduction of compensating networks. 

On the completion of the initial tests the equipment was applied to the 
measurement of the Cygnus A and Cassiopeia A radio sources. A typical record 
is shown in Fig. 4. The different declinations of the two sources necessitated 
the use of slightly different settings of the artificial delay lines and rotating phase 
shifters. At most aerial spacings it was possible to use a setting for the rotary 
phase shifters such that both sources gave fringes of suitable period but of opposite 
angular velocity. The direction of the velocity vector determined the interpre- 
tation of an observed phase shift as leading or lagging, and hence the sense of the 
asymmetry of the source. 
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FIG. 4. 


It was found that the measurement of the argument of the fringes was a 
considerably more powerful method for the determination of the order of 
the maxima of the transform patterns than a similar small number of measurements 
of the fringe modulus. The positions of the minima of the Cassiopeia source 
could be most accurately assessed by observing the change of phase in their 
immediate vicinity. 

The measurements made on the Cygnus source confirmed the bifurcation 
of this object (I.A.U. Symposium on Radio Astronomy, 1955), whilst the 
measurements on the Cassiopeia source established three maxima of the Fourier 
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transform of its intensity distribution and showed the existence of a faint 
component of larger angular diameter off-set from the principal component 
of the source. The details of this work will be described in subsequent papers. 
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PESIOME JIOKJIAJIOB, B ITIEPEBOJIE HA PYCCKHM A3bIK 


MHTEHCHBHOCTH, TOJIAPH3SAUMA WU SJIEKTPOHHAA ITJIOTHOCTb 
COJIHEMHOM KOPOHBI TIO POTOCHUMKAM IIOJIYYEHHbIM IIPH 
TIOJIHOM COJIHEMHOM 3ATMEHHM 25-rO @PEBPAJIA 1952 I 


X. gon Karobep 


Kopona dbotorpadupoBasach Mp MOUHOM COJIHEUHOM 3aTMCHHH 25-ro despasia 1952 r. c 
yueTBepeHHOH Kamepoili, c (OKYCHbIM paccTOAHHeM B 200 cM. [lomyueHo mIeCcTb CHAMKOB, Ka?K Abii 
C UeTbIPbMA H300paxKeHHAMH KOPOHBI, OHOBPCMECHHO IKCIOHHPOBaHHBIX, Tpoe Yepe3s NOJAPOMAbI 
C yIJIOM Me@?KJ{y MJIOCKOCTAMH NOApHsayHH B 60°, a ueTBepTbIii Ges NOMOIIM NoAponsa. 
Kaok bili CHHMOK CHaG>KeH OTHOCHTeJIbHOM HM aGcouOTHOM KasMOpaynamu. IIInpokoe mpHmeHe- 
HHe (OTOMeTPH4eCKHX H3MeEPeHH NOSBOJIMIIO OIpeeJIMTh HHTCHCHBHOCTS, BeJIMUHHY HM HalipaB- 
jieHHe MOMApHsaljMH, OTAeIMTA KOpoHbI K oT F' u omnpeyesIHTb 3JICKTPOHHbIC MJIOTHOCTH, 
PesysbTaTbI PeSIOMHPOBaHb! rpaduueckH. OHM cpaBHHBalOTCA C MOJeJIbIO BaH fle XIOUBCTa HC 
PyTHMu HaOJWOeHHAMH TOFO 2Ke 3aTMeHHA. 


SIKBUBAJIEHTHBIE KBAHTOBbIE CBETOBbIE OTIJAUM POTOSMYJIbCHH 
Ilémp Deavzeem 


BbIBOATCA SKBHBAJICHTHOe KOJIMUECTBO HAKOMICHHbIX (POTOHOB H 9KBHBAJICHTHbIC KBAHTOBbIe 
CBeTOBbIe OTWAUH JIA “eTLIPex rpynN oTosMyJIBCHH H3 OAHOPOAHOrO cocTaBa MO YyBCTBHTeJIb- 
HOCTH HM MO 3epHuHcTocTH. Bsicume oOHapyxKeHHbIe MOHOXpoMaTHUeCKHe OTQZaWH=0,9% 
TIpH CHHeM IyBeTe (430 Mu) 4 0,6% mpH KpacHom (607 Mu). Bpicuime oTWaun OKaSbIBaIOTCA TIpH 
IIOTHOCTH TIpeBbiuaroujeH BCcero Ha O,I WO 0,2 MIOTHOCTh ByamH. Bpicuiee 9KBHBaIeHTHOe 
KOJIMYECTBO HaKOMJICHHbIX (POTOHOB ABJIAeTCA B Mpeyeuax OT 0,03 Ha uw? Aa “ ObicTpeimen ” 
QMYJIBCHH MpH CHHeM UBeTe (430 Mu) WO 0,4 Ha mu? JIA CaMOTO MeJIKO3eEPHHCTOrO MaTepHasia 
pH KpacHomM Bete (607 My). TloquepkupaeTcH paSHHlla M@KIY BbIqepoKKaMH A 
MaKCHM@JIBHOM OTWauH HW JIA MAKCHM@JIBHOrO 9KBHBAJICHTHOTO HaKOMIeHHA. 


BEPOATHOCTH MEPEXOJIA 3AIIPEMEHHbIX JIMHHM Nin UW Nim 
P. X. Tapcmane 


BbIunCHAIOTCA BEPOATHOCTH Mepexofa WIA MSIIy4eHHii MarHUTHBIX AuMONeH HW 9IeKTpH4UeCKHX 
KBaypynouei WIA BCex MepexO0B M@KZy yPOBHAMH KOHdurypayun 3d° noHoB Nilll 4 Me*KLy 
YPOBHAMH KoHdurypaynit 3d*4s u 3d° nonoB Nin. 


BbIMUUCJIEHHE PAKTOPOB TAYHTA JIJIA CBOBOJHO-CBOBOJIHbIX 
MIEPEXOJIOB BBJIMNSH TOJIOXKMTEJIBHbBIX HMOHOB 


H. ®. Fpaum 


OGcyxWaloTCA MeTOAbI MOryliwe CiIyKHTb JIA BbIUHCIeHHA dakTOpoB [Taynuta AA 
CBOGOHO-CBOGOZHBIX MepexoqoB BOmM3sH sm0O0ro MomOKMTeIbHOrTO HOoHa. JoKasbiBaeTca 
OKBUBaJIeHTHOCTh dbopmys Commepdenboa u Menyena u Ilekepuca aa BogOpoqonogq0bHBIx 
CHCTeM, a TakoKe OOCYKTaIOTCH HECKOJIBKO TIpHOMKeHHI K pelieHu1o sTHx cbopmyn. JioKasbi- 
BaeTCAH, UTO MO?KHO ObICTPO OMpeeIMTb UMCIICHHbIeC BeJIMYMHbI MOHOXPOMATHUeCKHX (paKTOPOB 
TayHta mpH s1060i CKOpOcTH 9JIeKTpOHa, C MOMOLIbIO OOLUINPHOM TaOJIMUbI BeIMYHH HM MpOcTBIx 
ACMMITOTHYUCCKHX pa3JIOxKeHHil. 

BumuanHe 9KpaHHpoBaHHA Ha cbakTopbi [ayHta HCUMCIIAeTCA JIA JJIMHHbIX BOJIH, IPHMeCHAA 
mpoctyro dopmysly, BbIBeeHHy!O 43 NpHOmmMKeHHA BopHa MpH NOMOWM yHKUMH aTOMHOTO 
nmoTeHuvasia B Bue r-TLonexp(—fnr). Mertog moscnaetca mpHmMepaMH, B KOTOPbIX yHKUHH 
Tomac-@Depmu HM CaMOCOrsmacOBaHHble NOTeCHIHaJIbHbIe (PYHKUMH BbIpaxKeHbI TaKHM OOpa3som. 

I[aetca mpocroii MeTO, BbIUMCIeHHA CpeqHHX BeMUNH dbakTOpoB TayHTa, B COOTBeTCTBHH 
paciipefeseHua ckopocteii no MakcBessly, NPHMeHAA KBaypaTypHy!o (OpMyJIy OCHOBaHHYIO Ha 
Ole€HKe MOAMHTerpasIbHOrO BbIpaxKeHHA B HMCXOJHBIX TOUKaX NOJMHOMOB JIareppa. TpxH TouKH 
MHTerpalMH OKaSbIBalOTCA JOCTAaTOUHbIMH JIA MOUyYeHHA TOUHOCTH pellieHHA CpaBHHMOH C¢ 
TOUHOCTBIO JOCTHTHYyTOM Gosee CO*KHbIMH MeTOaMH, 





Russian summaries of preceding papers Vol. 118 


ATOMHBIE HW MOJIEKYJIAPHbIE OTPHUATEJIBHBIE HMOHBI B 
ATMOC®@EPAX 3BE3]I 


JI. M. Bpancrom6 u B. 3. To. [ladonea 


B ciiekTpax 3Be3, NO3HHX KJIaccOB HaOOaeTCA HeMpepbIBHOe MOTIOUIeHHe, He CBASAHHOE C 
voHamMw H-. Wcrounnku 9Toro NOrsiouleHHA BMepBbie NOABIIAIOTCA Ha CosIHye, a MX MHTCHCHB- 
HOCTb YBCJIMUMBaeTCA MIpH Mepexoye K 3Be3qam Goslee MO3MHHX CIICKTPOBbIX KJIACCOB, B 
oco6eHHOCTH 3Be34aM yruiepoqHOoro Tua. OcHOBbIBaACh Ha H€aBHO MOJYYeHHbIX TeOpeTH4UeCKHX 
HM JladopaTOpHbIX JaHHbIX, KaCalOuIMXCAH aTOMHbIX HM MOJICKYJIAPHBbIX OTPHIaTeIbHbIX HOHOB, 
KOHIeHTpauun C~, O-, OH- nu CN- knaccudbuumpyrotca 06MM oO6pa3so0M Kak rHvraHTbi HM KapsIMKH 
Pa3JIM4HOrO COCTaBa HM OMpeesineTCA BIIMAHHe HX Ha KOI*PHUKeHT HeMpepbIBHOTO NOTIOUICHHA. 
O- Kak paBuHyio HesHauMTembHa; C- BepOATHO MMee€T BayKHOe 3SHAYeHHe y 3Be3j{ C OUeCHb 
OOJIbIIHM COJep»KaHHeM yruiepofa, KaK HampHmMep R Cor Bor; OH- mookeT uMeTb 3Ha¥eHHe y 
3Be3, NO3QHHMx KyIaccop M u S; CN~- BosmorKHO Gosee BaxkHa 4uem Cg TPH NOroujeHHH B 
@uoueTOBOM HallasoHe y 3Be3 mosqHero Kyacca N. Jina mposepku HaOsOgeHHeM 9OTMX 
Pe3yJIbTaTOB HEOOXOAHMbI TOUHbIC ClieKTPO*OTOMeTPHUeCKHe H3MePeHHA B WIMPOKOM JMana3soHe 
AJIMH BOJIH HW [ajibHeiaa paOota B JaGopaTopHH. 


O MIPEOBPASOBAHHH H3MEPEHHM MOJIAPH3AUMH MCKAXKEHHbIX 
IIOJIAPHSALIMEN TIPHBOPOB 


A. Jon. Becceaunx 


PaccMoTpeHO BJIMAHHe HaKJIOHHOrO OTpaxKeHHA Ha M3MepeHHA NONApHsalHH. BpiBosATcA 
Mmpoctpie PopMyJIbI, GuarosapxA KOTOPbIM MO?KHO HCMOJb30BaTb HaOJOAeHHA Haj 3Be3jaMH, Ubu 
NOJAPH3allMH M3BeCTHbI, WWI KaliMOpHpoBaHHA HaOmOTeHH NOABep>KeCHHbIX HaKJIOHHOMY 
oTpaxeHuio. Meroy o0bacHsAeTA HaOsOReCHHAMH r->KH E. pan II. Cmut. O6painaetca BHuMaHHe 
Ha IIPHTOJ{HOCTb BCJIMUHH Pp cos 26, p sin 20 npu pa3s00pe BHYTpeHHeM TOUHOCTH HM TIpH comocTaB- 
JIGHHM He€3aBHCHMbIX PAOB. 


MCNOJIBZOBAHHE HHTEP®EPOMETPA UYBCTBHTEJIBHOTO K 
M3MEPEHHAM @A3bI JIA USMEPEHHM BbIPAXKEHHH ®YPbE 
TIPOCTPAHCTBEHHBIX PACIIPETEJIEHHMM ACHOCTM TPH MAJIbIX 
YIJIOBbIX CMEIJEHMVAX 


P. C. [[owennucon 


OnucbipaeTcA MeTON, CyxKallMii QA OfMHAapHOrO oMmpeyeueHHA aMIVJIMTyQbI UW (pasbi 
KOMIUICKCHOrO BbiIpaxkKeHHA Cypbe yUHTbIBaIOWero pacpefeseHHA ACHOCTH B MpoOcTpaHcTBe 
IPH MaJIOM yIJIOBOM CMeLIICHHH, MOJIyYeHHOrO MOMOLIbIO pxAZla H3MepeHHi TPpoMHOM CHCTeMbI 
HuHTepdbepometposB. Lismepsetca aOcounOTHaAH BeJIMUMHa aMIVJIMTyQHOM dyHKUMH, Tora Kak 
casa M3MePHeETCH OTHOCHTCJIbBHO JaHHOM BeJIMUMHbI, NOyYaemOH MpH KOPOTKOM paccTOAHHH 
M@*Ky aHTeHHaMH. ONIMCbIBaeTCA OMbITHbIM pagHOuacTOTHbIa HHTepdepomeTp, MCIMOJIb3SyIOUMH 
9TOT NPHHUKN HW pasOupaeTtca ero ZelicTBHe. 
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